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FOREWORD 

Work  described  in  this  report  was  accomplished  by  the  Flight  Mechanics 

Division  of  the  Air  Force  Flight  Dynamics  Laboratory  and  the  Digital 

programming  section,  4950th  Tesc  Wing  under  Project  1431,  "Flight  Path 
Analysis,"  Task  143109,  "trajectory  and  Motion  Analysis  of  Tlight 
Vehicles."  The  formulation  and  interim  documentation  were  completed 
by  Major  Urban  h.  D.  Lynch.  Programming  was  accomplished  by  Mr.  Fay  0. 
Young  of  the  Digital  Programming  Section  (ADDP),  Computer  Science  Center 
4950th  Test  Wing. 

This  report  was  prepared  by  Mr.  John  J.  Dueweke  of  the  High  Speed 
Aero  Performance  Branch  (FXG),  and  Mr.  Fay  0.  Young,  and  combines  the 
l applicable  portions  of  FDL-TDR-64-1 , Part  I,  Volume  1,  with  vhe  interim 

[ documentation.  The  overall  report  is  divided  into  four  parts: 

Part  I.  Capabilities  of  the  Takeoff  and  Landing  Analysis  Computer 
Program 

Part  II.  Problem  Formulation 
Part  III.  User's  Manual 
Part  IV.  Programmer's  Manual 

This  report  was  submitted  by  the  authors  in  March  1972. 
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ABSTRACT 

A well-defined  integration  of  the  various  aspects  of  the  aircraft 
takeoff  a’nd  landing  problem  is  presented  in  the  form  of  a generalized 
computer  program.  Total  aircraft  system  performance  is  evaluated  during 
the  glide  slope,  flare,  landing  roll  and  takeoff. 

The  flight  dynamics  of  a generalized,  rigid-body,  aerospace  vehicle 
are  formulated  in  six  degrees  of  freedom.  A flat,  nonrotating  Earth  is 
assumed.  The  independent  equations  of  motion  of  up  to  five  oleo-type 
landing  gears  are  also  formulated. 

A control  management  formulation  is  developed  to  automatically  adjust 
control  variables  to  correct  errors  in  the  "ehicle's  dynamic  state. 
Stability  in  the  small  is  used  to  maintain  stability  in  the  large. 

The  equations  of  motion  are  integrated  using  a generalized  variable- 
step  Runge-Kutta  technique. 

The  formulation  is  programmed  for  the  CDC  CYBER  70  and  6000  series 
computers  using  the  SCOPE  3.4  operating  system.  The  entire  program 
is  written  in  Fortran  Extended 
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SECTION  I 
INTRODUCTION 

In  the  design  of  an  aircraft,  the  engineer  is  confronted  with  the 
problem  of  takeoff  and  landing  and  the  design  of  aircraft  systems  and 
techniques  to  perform  this  function.  The  final  evaluation  of  these 
systems  lies  in  the  answer  to  the  question:  How  does  the  aircraft  and 

its  systems  perform  as  a unit?  The  Takeoff  and  Landing  Analysis  (TOLA) 
Computer  Program  is  the  result  of  an  attempt  to  generalize  the  aircraft, 
the  main  aircraft  control  systems,  and  the  landing-takeoff  situation 
into  a single  comprehensive  calculation  to  answer  this  question. 

The  TOLA  simulation  answers  the  above  question  in  the  form  of  a 
well-defined  integration  of  the  various  aspects  of  takeoff  and  landing. 
In  the  equations  of  motion  the  assumption  is  made  that  the  main  aircraft 
frame  is  rigid;  however,  the  dynamic  effects  of  up  to  five  independent 
landing  gears  are  included  in  the  equations.  The  position  and  velocity 
of  each  strut  and  secondary  piston  are  obtained  by  numerical  integration 
subject  to  position  constraints  (for  example,  the  main  strut  must  move 
within  the  limits  of  the  fully  extended  position  and  strut  bottoming 
position).  The  same  form  of  solution  applies  to  the  aircraft  itself. 
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SECTION  II 

COMPUTER  AND  SYSTEM  REQUIREMENTS 

The  Takeoff  and  Landing  Analysis  Computer  Program  (TOLA)  has  been 
written  for  use  with  the  CDC  CYBER  70  and  6000  series  computers  using  the 
SCOPE  3.4  operating  system.  The  entire  TOLA  is  written  in  Fortran 
Extended . 

1.  CYBER  70  or  6000  series  computer 

a.  A CYBER  70  or  6000  series  computer  with  32K  (decimal),  or 
larger  core. 

b.  Six  CDC  tape  transports. 

c.  Control  Data  Card  Reader. 

d.  Control  Data  Printer. 

e.  12  inch  CALCOMP  off-line  plotter. 

2.  The  CDC  tape  transports  may  be  replaced  by  other  equipment  which 
will  simulate  magnetic  tapes  such  as  disk  storage,  except  for  one  tape 
unit  that  may  be  used  to  generate  a plot  tape  for  the  CALCOMP  off-line 
plotter. 
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SECTION  III 
PROGRAMMING  CONCEPTS 

1.  THE  USE  OF  COMMON 

Whenever  possible,  a variable  is  placed  in  the  FORTRAN  "COMMON"  Area. 
There  are  several  reasons  for  this: 

a.  The  communications  between  subroutines  is  simplified. 

b.  The  structure  of  the  directory  is  simplified.  Since  the  number 
of  variables  i r,  COMMON  is  quite  large,  all  COMMON  cards  are  not 
placed  in  each  assembly/compilation.  Instead,  required  "duirmy" 
cards  are  placed  in  each  deck  of  source  cards.  This  has  in  a 
small  manner  reduced  the  number  of  COMMON  cards  in  each  deck. 

2.  TABLES  AND  TABLE  USAGE 

One  of  the  usual  required  modifications  of  any  program  is  the  change 
of  table  sizes.  With  this  in  mind,  a COMMON  block  of  locations  has  been 
set  aside  and  the  required  number  of  ceils  for  each  tabie  is  specified 
with  data  (see  TABRE  for  data  preparation).  This  requires  no  reassembly 
or  recompilation  unless  the  total  number  of  cells  required  exceeds  tne 
COMMON  block  of  800  cells. 

3.  SYMBOLIC  INPUT 

Although  the  FORTRAN  system  itself  has  a system  of  input  routines, 
the  program  does  the  actual  translation  of  the  cards  using  special  coded 
routines.  Input  data  may  be  read  using  a system  of  symbols  which  is 
designed  to  give  engineering  meaning  to  the  analyst.  The  symbols  are 
referenced  to  actual  locations  by  the  use  of  COhWON  and  subscripts. 
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4.  TRAJECTORY  PRINTING  METHOD 

The  printing  of  a trajectory  may  be  divided  into  four  categories, 
a.  Initial  Printing  - The  printing  of  specific  values  at  the 
first  stage  and  at  each  subsequent  major  stage. 

> b.  Code  Printing  - The  printing  of  codes  which  will  identify 

| the  variables  which  are  to  be  obtained  in  the  coming  time  history  print. 

[ c.  Time  History  Printing  - The  printing  of  values  specified  at 

the  requested  points  of  the  trajectory. 

d.  Diagnostic  Error  Printing  - The  printing  of  errors  detected 
by  the  program. 

All  input  data  involved  for  a case  is  printed  on  the  output  page 
preceding  the  computation  of  the  first  stage  printout.  Also,  data  read 
in  at  stage  times  will  be  printed  out  between  the  stages  of  the  trajectory 
output. 

Initial  print  is  designed  to  print  certain  values  which  will  be 
constant  during  the  trajectory  and  serves  as  a reminder  of  what  valuer 
have  been  used  for  these  constants. 

Code  printing  is  performed  once  per  major  stage  to  identify  the  time 
history. 

The  time  history  print  is  designed  to  print  in  a minimum  space. 

That  is,  if  a certain  variable  is  not  desired  as  output,  it  is  not  printed 
1 and  other  desired  variables  are  moved  in  the  print  format  accordingly. 
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The  entire  printing  is  controlled  to  print  on  a page  11  x 14  inches 
and  will  print  a maximum  of  51  lines  per  page.  Page  ejection  and  lines 
control  are  provided  by  the  subroutines  DEF  and  LINES. 


TAPE  USAGE 

Tapes 

Equipment 

Usage 

Tape  5 

Disk  or  Tape 

Data  Input 

Tape  6 

Disk  or  Tape 

Printed  Output 

Tape  13 

Disk  or  Tape 

Data  saved  to  be  used  to 

generate  a plot  tape  by 

Plot  program  (PLTSDF) 

Tape  16 

Disk  or  Tape  I 1 

| Used  by  the  symbolic 

Tape  31 

Disk  or  Tape  ( , 

| input  routine  to  save 

1 

i input  data 

Tape  7 

Tape 

CALC0MP  Plot  tape 

The  above  describes  the  tape  usage  other  than  for  the  FORTRAN  system. 
All  modification  of  tapes  required  may  be  made  with  control  cards  placed 
in  front  of  the  program  before  submitting  to  the  computer. 

6.  STRUCTURE  OF  PROGRAM 

Due  to  core  storage  limitations  (32K),  it  was  necessary  to  use  the 
Overlay  feature.  The  following  is  the  structure  of  the  program: 

a . OVERLAY  (0,0)  (FOR.  EXT) 


TOLA 

STFL 

TRNPOS 

FXE 

STFLD 

HI  HO 

INUPD 

S TOVAR 

TLU 

LN'JPD 

ARRAY 

TFFS1 
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INPUZ 

LINES 

IUTEG 

ASIN 

UPDAT 

ACOS 

VPCS1 

MIM1N 

A TAN  2 

SAC  SI 

lcdet 

ERROR 

aep.oi 

STGTS I 

LXERR 

AOL’AD 

STGTST 

NDTLU 

OPT1 

TDATA 

ATMS 

LG EAR 1 

INVR3 

LGEA3C 

ASRCU 

MULT 31 

SDFI.GP 

DEF 

OVERLAY 

(1,0)  (FOR  EXT) 

TOLAN1 

TSRCH 

TFFS2 

DSERCH 

READ 

PACKRR 

DIPLAC 

RITE 

TAhRE 

D I ROD A 

RI'.ADA 

DIR] DA 

STORE 

DIR2DA 

WRCARD 
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c.  OVERLAY  (2,  0)  (FOR  EXT) 


TOLAN2 

ENGREV 

AUTS 

ENGFL 

FLARE 1 

CENCL 

AUTPR1 

CTENGL 

THAUTS 

TFFS8 

TFFS9 

7.  PROGRAM  ORGANIZATION 

The  TOLA  Computer  Program  is  written  in  FORTRAN  Extended.  The 
program  is  segmented  and  takes  advantage  of  the  FORTRAN  overlay  features. 

This  section  attempts  to  describe  the  overall  organization  of  the 
program  from  the  viewpoints  of  control  cards,  tape  usage,  deck  set-up, 
and  organization. 

The  program  is  broken  up  into  three  overlays  as  follows: 

a.  Overlay  (0,  0).  Contains  all  system  routines,  main,  executive, 
integration,  computation  of  the  equations  of  motion,  and  printing. 
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b.  Overlay  (1,  0).  Set  up  tables,  input  routines,  and 
input  directory. 

c.  Overlay  (2,  0).  All  routines  of  the  Autopilot. 

Plotting  tapes  are  generated  by  a separate  program  for  plotting  on 
the  Cal  Comp  plotter. 

a.  Storage  Reference.  All  variables  requiring  arrays  have  been 
arranged  in  the  standard  FORTRAN  convention;  for  example,  an  array  A^  is 
stored  in  increasing  storage  locations  for  increasing  i.  Matrices  are 
stored  columnwise. 

b.  Integers.  All  integers  are  assumed  to  be  in  a 60-bit  word, 
right-justified. 

c.  C^MM^N . In  order  to  decrease  the  length  and  time  required  in 
calling  sequences,  liberal  use  of  labeled  C0MM0N  has  been  made.  For  the 
actual  variable  and  their  arrangement  in  C0MM0N , the  user  is  referred 
to  the  program  listing. 

d.  Variable  Names.  Because  any  variable  may  be  referred  to  by 
FORTRAN,  all  integer  variable  names  begin  with  the  leading  letters  I,  J, 
K,  L,  M,  or  N.  This  does  not  mean  that  all  noninteger  variable  names 
begin  with  letters  other  than  I,  J,  K,  L,  M,  or  N.  They  may,  in  some 
subprogram,  be  declared  integer  or  real. 

8.  DECK  SETUP 

a.  Running  the  TOLA  Program  requires  a particular  deck  setup.  The 
deck  structure  is  presented  as  a guide  only  in  determining  this  se^up. 
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b.  CONTROL  CARDS  (CDC  CYBER  70  or  6000  Series,  SCOPE  3.4).  All 
control  cards  are  left  justified  in  column  1.  The  end  of  record  is  a 
7,  8,  9 punched  in  column  1 and  on  end  of  job  card  is  a 6,  7,  8,  9 punch 
in  column  1.  In  the  control  card  examples  an  e.id  of  record  and  an  end  of 
job  will  be  used  in  place  of  the  cards. 

(1)  The  following  control  cards  will  execute  the  TOLA  Computer 
Program  from  an  UPDATE  tape  and  not  print  a listing  of  TOLA. 

Job  Card 

LABEL,  0LDPL,  R,  L = T0LACP , VSN  = tape  No.  RING  0UT 
UPDATE,  F. 

FTN,  I = COMPILE,  L = 0. 

RETURN,  0LDPL. 

LDSET,  PRESET  = ZERO. 

L0AD,  LG0. 

N9G0. 

TOLA. 

End  of  Record 

Changes  to  T9LACP  if  any  in  UPDATE  format 
End  of  Record 
Data  Cards 
End  of  job. 

If  a listing  Is  desired,  omit  the  parameter  L = 0 on  the  FTN  control 
card. 
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(2)  The  following  control  cards  will  generate  a new  UPDATE  tape,  an 
absolute  file  on  tape,  and  list  the  T0LA  Program. 

Job  Card 

LABEL,  0LDPL,  R,  L *»  T0LACP,  VSN  = tape  No.  RING  0UT 
LABEL,  NEWPL,  W,  L = T0LACP,  VSN  = tape  No.  RING  IN 
UPDATE,  N,  F. 

FTN,  I = COMPILE. 

RETURN,  0LDPL. 

UNLOAD,  NEWPL, 

LABEL,  TOLA,  W,  L * T0LACPABS,  VSN  = tape  No.  RING  IN 
LDSET,  PRESET  = ZERO. 

LOAD,  LGO. 

NOGO. 

TOLA. 

End  of  record 

Changes  to  TOLA  if  any  in  UPDATE  format 
End  of  Record 
Data  cards 
End  of  job 

(3)  The  following  control  cards  will  execute  from  an  absolute  file 
on  tape,  TOLA  is  an  absolute  file  on  a tape. 

Job  card 

LABEL,  TOLA,  R,  L = TOLACPABS,  VSN  = tape  No.  RING  OUT 
TOLA. 

End  of  record 
Data  cards 
End  of  job 
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(4)  If  it  is  desired  to  save  data  on  tape  for  CALC9MP  plotting, 
include  the  following  control  card  with  the  LABEL  cards: 

LABEL,  TAPE  13,  W,  L = T0LADATA,  VSN  = tape  No.  RING  IN 

(5)  The  following  control  cards  will  generate  a new  updated  program 
on  permanent  file  ( PF) , generate  on  absolute  file  of  T9LA  on  PF,  and 
execute  T0LA. 

Job  Card 

LABEL,  9LDPL,  R,  L = T0LACP,  VSN  = tape  No.  RING  OUT 
REQUEST,  NEWPL, *PF. 

UPDATE,  N,  F. 

RETURN,  9LDPL. 

CATALOG,  NEWPL,  T0LACP,  RP  = 999,  CY  = 1,  ID  = Prob  No. 

RETURN,  NEWPL. 

F TN,  I = COMPILE,  L = 0 
REQUEST,  T0LAP,  *PF . 

LDSET,  PRESET  = ZERO. 

LOAD,  LGO. 

NOGO. 

TOLA. 

CATALOG,  TOLAP,  TOLACP,  RP  = 999,  CY  = 2,  ID  = Prob  No. 

End  of  record 

Changes  to  TOLA  if  any  in  UPDATE  format 
End  of  record 
Data  cards 
End  of  job 
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(6)  To  execute  from  a permanent  file,  use  the  following  control 
cards: 

Job  Card 

ATTACH,  TOLA,  T0LACP,  CY  = 2,  ID  = XXXXXX. 

TOLA. 

End  of  record 
DATA  cards 
End  of  job 

(7)  To  execute  from  a permanent  file,  and  generate  a plot  tape 

for  the  CALC0MP  plotter  on  TAPE  7,  use  the  following  control  cards.  The 
plot  tape  generating  program  (PLTSDF)  is  located  on  PF:  PLTSDF,  CY  = 1 . 

Job  Card 

ATTACH,  TOLA,  TOLACP,  CY  = 2,  ID  = XXXXXX. 

TOLA. 

RETURN,,  TOLA 

REQUEST,  TAPE  7,  MT,  HI,  N,  VSN  = Tape  No.  Ring  IN 
ATTACH,  PLTSDF,  PLTSDF,  CY  = 1 , ID  = XXXXXX. 

PLTSDF. 

End  of  record 
Data  cards  for  TOLA 
End  of  record 
Data  cards  for  PLTSDF 
End  of  job 
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9.  DATA  FORMAT 

Card  Format  - The  program  input  routine  (READ)  expects  the  following 
format. 

. _ I II  a ■ • t 

73  - 80 
VIII 

Card  Field  I - Contains  the  symbolic  name  of  the  variable  which  data 


Card -Columns  - 

1 - 6 

7 

8-10 

11  1 12  - 66 

67  - 72 

Field 

I 

II 

III 

IV  1 V 

VI 

contained  in  Field  V begins  loading.  Example:  Card  Column 


1 

IGAM7D 

SIG7D 


12 

-1.23 

90. 
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Card  Field  II  Not  Used 

Card  Field  III  - Contains  the  words  DEC,  0CT,  BCD,  TRA,  INT,  or  is 
blank  depending  on  the  type  of  data  to  be  loaded.  The  word  OCT  indicates 
that  the  data  is  to  be  interpreted  as  octal  numbers.  The  word  BCD 
specifies  that  N binary  coded  decimal  words  (N  punched  in  column  12) 
beginning  in  column  13  are  to  be  loaded.  The  word  TRA  denotes  to  the 
input  routine  that  all  data  has  been  read  and  to  return  control  to  the 
calling  program.  The  word  DEC  and  blank  are  equivalent  and  specify 
that  data  loaded  is  decimal  data. 


0CT  Example 
Card  Column 

1 

8 

12 

NSMAIN 

OCT 

17 

BCD  Example 
Card  Column 

1 

8 

12 

REM 

BCD 

3SDF2-GEAR-M0D 

The  3 in  Column  12  specifies  3 words  where  each  word  is  considered  to 
be  6 characters  including  blanks.  The  largest  number  of  6-character 
words  that  can  be  loaded  from  one  card  is  9.  The  analysts  should  be 
very  careful  to  see  that  the  BCD  information  does  not  get  punched  into 
Field  VI.  This  will  cause  an  input  error. 

DEC  Example 

Card  Column  1,8  12 

VTAB01  DEC  2,0., 1.67, 20000., 1.67 

Note  that  the  first  character  in  Column  12  is  an  integer;  the  input 

routine  will  load  only  one  integer  per  DEC  card,  and  that  has  to  be  the 

first  number  punched  in  Field  V. 

VTABOI  DEC  2., 0., 1.67, 20000., 1.67 
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If  the  above  card  is  punched,  the  two  will  now  be  loaded  into  the  machine 
as  a binary  floating  point  number.  The  other  numbers  will  be  loaded 
the  same,  with  the  decimal  point  assumed  right-justified. 

If  anything  other  than  0CT,  BCD,  INT,  TRA,  or  blank  appears  in 
Field  II  then  the  word  DEC  is  assumed. 

INT  Example 


Card  Column 

1 

8 

12 

IP 

INT 

1 

IP 

1, 

IP 

INT 

1, U,1 

When  the  word  INT  is  useci  it  is  assumed  that  all  numbers  on  the  card 
will  be  loaded  as  integers.  If  only  one  integer  is  punched  per  card  the 
INT  may  be  punched  or  omitted. 

Card  Field  IV-Not  Used 
Card  Meld  V 

The  actual  input  data  to  the  program  is  punched  in  the  Field  V. 

DEC,  INT  and  0CT  numbers  must  always  be  left-adjusted;  that  is,  it 
must  always  start  in  column  12  on  the  input  card.  All  numbers  are 
separated  by  a "comma"  and  the  field  terminates  with  the  first  blank. 

BCD  information  begins  in  Column  13  and  the  maximum  number  of  6-character 
words  per  card  is  nine.  Note  that  since  Field  V ends  with  the  first 
blank,  the  user  may  punch  any  comments  in  the  remainder  of  the  field. 

Card  Field  VI 

This  field  specifies  the  initial  subscript  of  the  data  in  Field  V. 

If  this  field  is  blank,  an  initial  subscript  of  1 is  Implied.  The 
subscript  may  appear  anywhere  within  the  field. 
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Example 
Card  Column 


1 

12 

PZER0 

30470.4,41538.24,41538.24 

PZER0 

42538.24,42538.24 

(or  blank) 


In  the  example  above,  the  number  30470.4  is  loaded  into  the  first  cell 
of  the  array  PZER0.  On  the  second  card,  42538.24  is  loaded  into  the  fourth 
cell  of  the  array.  The  one  and  four  punched  in  Field  VI  indicate  the 
subscript  for  the  array  P2.ER0. 


Card  Field  VII 

Not  used  as  far  as  the  input  routine  is  concerned.  This  may  be  used 
as  a sequence  number  for  the  card. 


10.  TABLE  FORMAT 

The  various  types  of  tables  used  by  the  program  may  be  classed  as 
follows: 


One  Dimensional  Tables 

Example  1:  NTIRE5 

Card  Column  1 

NSTRUT 

NTIRES 


n.  • f (i),  i = I,  2,...,  NSTRUT 

12 

5 

14. ,5. ,6. .6. ,6. 


INSTRUT=  Fixed  point  number  which  is  the  number  of  struts  on  the  aircraft. 

For  example,  the  number  of  tires  on  strut  2 is  6}  i.e.,  - f ( 2 ) = 6. 

i = independent  variable  values 

n.j  = Corresponding  dependent  variable  values 

Example  2:  Aerodynamic  Data 

Card  Column  1 12 

INDA01  1 

ATAB01  .0065,. 00748 
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INDAOI  jM  designates  that  there  are  data  in  ATABOl.  The  first  data  point 
is  for  full  ground  effect;  the  second  data  point  is  for  no  ground  effect 
in  all  aerodynamic  tables. 


Two  Dimensional  Table 

Example: 

Card  Column 


VTAB01 

1 

VTAB01 


12 


xcc=  f(m) 


N,M1  *XC61,Mz*XCG2,‘**,MN,XCGn 

N = Fixed  point  number  equal  to  2 times  the  number  of  independent  variables, 
For  a 20-point  table,  N would  equal  40.  The  total  number  of  machine 
cells  required  for  this  table  is  41. 

M.j  = Independent  variable  values 

Xrr  = Corresponding  dependent  variable  values 


N - Dimensional  Table 


Example: 
Card  Column 


T 

1 


F(N,K  ) 
12  N 


IT10W 

NN 

IT10X 

NMN 

TTAB10 

N1,N2»N3,',,,NNN 

TTAB10 

MN1,MN2,MN3,*‘,,MNNMN 

TTAB10 

TNrMNl,TN2,MNl  ’•••»TNNN,MN1 

TTAB10 

TN1 ,MN2,TN2,MN2 TNNN*MN2 

TTAB10 


TN1  ,MNNMN’TN2’MNNMN*-  • ' ,TNNN,MNNMN 


NN  and  NMN  are  fixed  point  numbers  of  independent  variables. 


TN,  . \ 


Tn  , M^^  are  values  of  independent  variables.  The  table  subscripts 
NN 

would  apply  to  the  N-dimensional  table  as  well  as  the  two  dimensional . The 
total  number  of  machine  cells  required  for  an  N-dimensional  table  equal 
NN  * NMN  + NN  + NMN. 
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Examples: 

C = F(X,  Y) 

Machine  cells  required 
C = F(X,  Y,  Z) 

Machine  cells  required 


NX  = 2 = points  for  X 
NY  = 2 - points  for  Y 
2x2  + 2 + 2 = 8 cel  1 s 


NX  = 20  = points  for  X 
NY  = 10  = points  for  Y 
NZ  = 15  - points  for  Z 
20  X 10  X 15  + 20  + 10  + 15  = 3045  cells 
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SECTION  IV 

FORTRAN  EXTENDED  OVERLAY  (0,0) 

1.  TOLA  - MAIN  PROGRAM 

a.  Purpose  - Initializes  parameters  in  the  read  routine  through 

COMMON,  initializes  table  cata,  and  parameters  in  storage  routines,  and 
calls  EXE. 

b.  Usage  - Calls  the  executive  routine  (EXE)  for  each  case. 


2.  EXE  - EXECUTIVE  ROUTINE 

a.  Purpose  - To  zero  all  variables  that  may  be  read  in  as  input, 
initialize  subprograms,  and  set  nominal  values.  Read  input,  do  all  the 
proper  initialization,  set  up  tables  dimensions,  check  for  staging, 
printing,  etc. 

b.  Usage  - The  executive  routine  is  the  controlling  program.  When 
a case  is  completed,  a return  is  made  to  TOLA. 


3.  STGTSI , STGTST  - STAGE  TESTING  ROUTINES 

a.  Purpose  - To  test  the  possiblity  of  staging  on  any  of  up  to  four 
increasing  and  four  decreasing  variables. 

b.  Method  - Given  the  four  increasing  variable  BCD  names  (in  arrcy 
AINCRS)  and  the  four  decreasing  variable  BCD  names  (in  array  DECRES)  the 
routine  first  searches  the  directory  for  their  location  (routine  STGTSI). 
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FLOW  DIAGRAM  — MAIN  P*’:OGRAM  (TOLA) 


AFFDL-TR-71-155 
PARI  IV 
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i 

! 
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FLOW  DIAGRAM  - EXECUTIVE  ROUTINE  (EXE) 


TPO  * T+DELTS 

TPO  < Tmo*  ) 

NO 

OELTS  * 
H0*  0 
END  * 

Tmo*‘T 

ELTS 

TRUE. 

PRINT  AUTOPILOT 
VARIABLES; 

0 1CONTR  « 2 
©CALL  OVERLAY  (2,0) 

#0 

INDAUT  ^ 


INTEGRATE  ALL  VARIABLES 
CALL  MIMIN  (MIM) 

(IF  MIM  #0,  PRIMIN 
HAS  BEEN  EXCEEDED.) 


0TIMEP  * T PRINT-1. E-6 
PRINT  SOF  8 
SUBPROGRAM  VARIABLES 
©CALL  0PT6 


SAVE  DATA  ON  TAPE: 
CALL  SDFL6P 


COMPUTE  AUTOPILOT 
VARIABLES  : 

0 ICONTR * 2 
©CALL  OVERLAY (2,0) 


1NDSTE 

I #0 


L±°_ 

INDAUT 

YES 
T <TIMEP 


TIMEA  ■ T + DELTS-I-  E-  6 
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FLOW  DIAGRAM  - EXECUTIVE  ROUTINE  (EXE.) 


1 
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In  routine  STGTST,  and  when  ISTAGE=0,  the  routine  may  then  be  used  to 
test  for  any  of  the  increasing  variables  being  greater  than  the  values 
in  the  STEST  array  (i.e.,  increasing  variable^  ^STESTi).  In  like  manner, 
the  decreasing  variablei  > STESTDi  test  is  made  for  1 = 1,  2,  3,  4. 

Note  that  the  testing  stops  at  the  first  test  to  be  satisfied.  If  a 
test  is  satisfied  for  increasing  variables,  the  routine  places  the  BCD 
name  of  the  varible  in  STGVAR  and  the  BCD  word  "INCR"  in  STGVAR-1.  This 
process  is  similar  for  decreasing  variables,  except  that  ,,DECR,,  is 
placed  in  STGVAR-1. 

Similarly,  when  ISTAGE  / 0,  the  routine  may  be  used  to  test  for  al 1 
of  the  increasing  variables  being  equal  to  the  values  in  the  STEST 
array  (i .e., increasing  variablei  = STESTi).  In  like  manner,  the  decreasing 
variablei  = STESTDi  test  is  made  for  i = 1,  2,  3,  4.  Note  that  in  order 
to  pass  the  test  the  conditions  must  be  met  for  all  i. 

c.  Usage 

(1)  Initialization  - The  statement,  CALL  STGTSI 

must  be  given.  This  statement  must  be  given  each  time  a new 
set  of  variable  names  is  to  be  tested. 

(2)  Testing  Staging  - To  test  if  staging  is  indicated, 

CALL  STGTST  ( INDSTG) 

The  routine  returns  INDSTG  = 0,  No  staging 

INDSTG  = 1,  staging  indicated 
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4.  INUPD,  LNUPD,  INPUZ,  INTEG,  UPDATE  - INTERFACE  ROUTINES  FOR 
INTEGRATION  ROUTINE 

a.  Purpose  - To  serve  as  an  interface  between  the  integration  routine 
proper  (MIMIN)  and  any  routine  requesting  a variable  to  be  integrated. 


There  are  five  logical  functions  which  these  routines  perform. 

For  a particular  call  to  one  of  these  routines,  one  of  these  functions 
will  be  enacted.  The  P array  is  the  array  of  current  derivative  values 
of  the  variables  that  are  being  integrated.  The  Y array  is  the  array 
of  the  current  integrated  variable  values. 


b.  Usage  - for  each  of  the  interface  routines. 

(1)  CALL  INUPD(N.L) 

The  number  of  integrated  variables  is  increased  by  N.  The  subscripts 
in  the  P and  Y arrays  for  the  values  XDOT  and  X respectively  are  stored 
in  the  array  L. 

(2)  CALL  LNUPD  (M) 

The  number  of  integrated  variables  is  decreased  by  M. 

(3)  CALL  INPUZ 

The  P an-1  Y arrays  are  set  to  0;  the  number  of  integrated  variables 
is  set  to  0. 
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(4)  CALL  INTEG  (K,XDOT) 

The  value  of  XDOT  is  stored  in  P(K). 

(5)  Call  UPDAT  (JX1,  JX2,  XJ1 , XJ2,  XJ3,  XJ4 , XJ5) 

JX1  = No.  of  variables  (1  <_  JX1  £ 5);  JX2  = subscript  of  first 
variable  (XJ1);  2nc*  variable  (XJ2)  has  subscript  of  JX2+1 , etc. 

When  the  logical  variable  3KIPUP  (DIRCOM)  is  false,  then  the  values 
of  the  integrated  variables  are  picked  up  from  the  proper  places  in  the 
Y array  and  put  in  XJ1 , XJ2,  etc. 

When  SKIPUP  is  true,  the  values  of  XJ1,  XJ2,  ....  XJ5  are  stored 
in  Y (JX2),  Y (JX2+1),  . . . , Y (JX2+JX1  - 1) 

INUPD  will  terminate  the  case  if  more  variables  are  requested  to  be 
integrated  than  there  is  room  for  in  the  integration  arrays;  at  the 
present,  this  upper  limit  is  50  variables. 

5.  MIMIN  - INTEGRATION  ROUTINE 

a.  Purpose  - To  perform  the  calculation  necessary  to  integrate 
an  arr$y  of  variables  by  the  variable-step  or  fixed  step  Runge-Kutta 
Method;  to  determine  an  estimate  of  the  relative  error,  and  from  this 
information,  determine  the  new  step  size  of  integration  in  the 
variable-step  mode  of  integration. 
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RETURN 


RETURN 
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b.  Usage  - The  entry  to  this  routine  is  as  follows: 

CALL  MIMIN  (SN) 

where  N = statement  number  to  which  nonstandard  return  is  made 
from  MIMIN. 

c.  Method  - The  variable  step  Runge-Kutta  Method  is  calculated  as 
follows: 

STEP  I 

Ymoxn  = | Yn  | 

Yn  = MX,Yn) 

YA«  = Yn  ♦ Yn 
YAn  s F (X  4 h/2  , VAn) 

VB„  = Yn  + “/2  YAn 

YBn  = F<X  + h/2  • W 
YCn  = Yn  ♦ h YBn 

YCn  2 F<X  + h«YCn> 

YPn  S Yn  ♦ h'6<Yn  + 2 YAn  4 2 f 4 Y-  ) 


II 

c 

>Q 

Yn  4 h/4  Yn 

YOn- 

F(X  + h/4,Y0n) 

II 

C 

Yn  + h/4  YD„ 

II 

c 

F (X  4 h/4  , YEfl) 

II 

c 

Y"  + h/2  YEn 

II 

F (X  4-  h/2  , Ypn  } 

Yln  = 

Yn  4 h/2  • 1/6  { Y, 

♦ YFJ 
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STEP  III 


YGn 

s 

Yin  + h/4  Yln 

YGn 

= 

F(X  + 3/4 h,  Y,  ) 
on 

YHn 

= 

Yln  + h/4  YGn 

YHn 

= 

F(X  + 3/4(it  YHn) 

YIn 

= 

Yin  ♦ YHn 

YIn 

= 

F (X  + h , YIn) 

Yn 

= 

Yln  + h/2  - 1/6  <Yln  + 2 Y6n  ♦ 2 YHn  ♦ Yj„) 

where  n = 1 , 2, 

3,  . . . , NDEFEQ 

NDEFEQ  = No.  of 

differential 

equations 

Computation  of  1 

Relative  Error 

ER  = ' 

iV|Y»  - W 

Zn  = Max 

(*-.„■  M 

) ; R = Mon 

[Rt  Hr  1/ Mo*  ( 

zn.  1 )] 

d.  Input  for  Integration 

Routine 

Symbol  Used 

Math 

Symbol  used 

Nominal 

by  READ  Routine 

Notation 

by  Integration  Value 

Remarks 

IVARBH 

IVARBH 

0 

Use  variable-step 

=1 , Use  Fixed  Step 

TIME 

t 

X 

0. 

Time  to  begin  integ. 

DELTS 

At 

ox 

.1 

Time  internal  to  int. 

AMINER 

At  mi  n 

DXMIN 

.001 

Minimum  At 

AMAXER 

At  max 

DXMAX 

10000. 

Maximum  At 

RELER1 

^EL 

tLFRl 

RELER1 

.00007 

Rel . error  tol.  #1 

RELER2 

rel 

tLER2 

RELER2 

.000005 

Rel . error  tol.  #2 

N 

N 

No.  of  diff.  eqs. 

Yi 

Y(N) 

Array  of  dep.  -var. 

Pi 

P(N) 

Array  of 

PRTMIN 

PRTMIN 

Print  Minimum 

INDLG 

Landing  gear  indicator 

INDSTE 

1 

ST0PIT 

SW 
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FLOW  DIAGRAM  (MIMIN) 


FLOW  DIAGRAM  (MIMIN) 
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6.  LGDET  - ROUTINE  TO  RESTRICT  LG  VARIABLES  IN  INTEGRATION  ROUTINE 

a.  Purpose 

When  the  statement  CALL  INTEG  (K,  XD)  is  executed,  the  derivative 
XD  is  stored  in  an  array  P at  location  K of  that  array  (i.e.,  P ( K )).  The 
integration  is  not  processed  until  a complete  pass  has  been  made  through 
the  program  and  all  calls  to  INTEG  have  been  made. 

As  the  integration  is  performed  the  Integrated  variables  are  stored 
in  an  array  Y at  the  correspondin'  location  K as  its  derivative  (i.e.,Y  ( K) ) . 

When  the  statement  CALL  UPDAT  (N,K,X)  is  executed,  the  integrated 
variable  X Is  transmitted  from  the  Y array  (location  Y (K) ) to  the  location 
X.  The  N designates  the  number  of  variables  to  be  transmitted. 

Due  to  the  requirements  of  the  landing  gear  problem  some  variables 
that  are  integrated  are  restricted  to  certain  conditions.  Therefore, 
it  was  necessary  to  write  the  routine  LGDET  to  restrict  these  variables 
in  the  integration  routine.  As  stated  above  these  variables  are  stored 
in  the  P and  Y array  in  the  integration  routine. 


b.  Linkage  - CALL  LGDET 
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FLOW  DIAGRAM  (LGDET) 


7L6PET 
■ENTRY , 


INDLG = 2 
♦ NO 

rTtr- 


1 RETURN) 


^Y(j)  >(-  ES(i))J 

t NO 

WRITE: " - ES(U  EXCEEDED" 
PRINT  VALUE  OF  S, 


1 = 1+1 
k = 21  + NSTRUT 
j * LA(k) 
kk  = k - 1 
jj  = LA  (kk) 

I - kk  + 2 it  NSTRUT 

II  * LA (£ ) 

w = k + 2*  NSTRUT 
mm  = LA  (nt) 
nn  = LA ( i ) 

HT1  = MT 


TTIME  = Y (j  )/|p(j)  | 


(ttime 

> ht  y- 

1 

| NO 

HTI  = TTIME 

f ycj)  < esi i ) y — ■ 
y—fHO 

(y))  ) <(SB(I)-ES(I)))— 

' f NO 

— (y^<(ES(I)  + SBU))) 
i NO 

write:  "es(I)  exceeded’ 

PRINT  VALUE  OF  S, 


Y(|j)  < O. 

I NO 


Y(jj)  = 0. 
P(j)  = O. 


1 

jf  >0 

t . 

TTIME  = (SBII)-Y(i))/P(j) 

1 

P()i  ) < 0. 
j^NO- 

P(ii)  = 0. 
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FLOW  DIAGRAM  (LGDET) 
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7.  ASRCH,  TDATA  - DIRECTORY  SEARCH  ROUTINES. 

a.  Purpose  - To  provide  a BCD  word  look-up  from  a subscript. 

b.  Method  - Given  a subscript,  the  routine  will  search  the  directory 
for  the  BCD  word  corresponding  to  that  subscript.  If  the  subscripts  do 
not  compare  the  BCD  name  is  set  to  blank  and  return  to  the  Calling  Program 
is  made. 

c.  Usage  - Entry  is  made  to  the  routine  with  the  following  statement 
a.  CALL  ASRCH  (LOCI,  SYM1 ) 

where 

LOCI  = subscript  being  searched 

SYM1  = the  variable  name  into  which  the  routine  is  to  store 
the  corresponding  BCD  name. 

d.  A call  to  the  subroutine  TDATA  initializes  table  name  data 
which  is  used  by  subroutine  TABRE  (in  Overlay  (1,0)). 
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8.  DEF  - HEADING  AND  PAGE  EJECT  ROUTINE 

a.  Purpose  - To  provide  page  ejection  and  title  printing. 

b.  Method  - Initially  the  current  page  number  (NPAGE)  is  incremented 
by  1.  The  title  is  printed  and  returned  to  the  calling  program. 

c.  Usage  - Entry  is  made  via  the  statement 
CALL  DEF 

9.  STFL,  STFLD,  STQVAR,  ARRAY  - Storage  and  Printing 
of  Output  Routines 

a.  Purpose  - To  provide  a method  of  printing  output  names  of 

variables  and  thei*1  values  when  necessary.  Names  of  variables  or  values 
which  are  to  be  printed  are  not  actually  printed  by  the  routines  until 
at  least  eight  have  been  accumulated  by  a series  of  calling  sequences. 

b.  Method  - Each  time  a call  is  made,  names  of  variables  or  values 

are  saved  until  eight  are  stored.  At  this  time  they  are  printed.  This 
process  is  repeated  until  all  names  or  values  have  been  handled.  If 
less  than  8 remain,  they  are  saved  for  further  call  statements  or  until 
forced  to  be  printed. 

c.  Usage 

(1)  The  printing  of  Hollerith  Code 
CALL  STFL  (JOPT,  N,  ARG1 ) 

JOPT  = 0:  Force  any  possible  remaining  print 

JOPT  = 1:  See  (2)  below  for  printing  of  Values  of  Variables. 

JOPT  = 2:  Prints  N words  of  Hollerith  information  from  ARG1  (1), 

ARG1  (2),  . . .,  AP.G1  (N)  where  there  is  at  most  6 characters  per  word. 
JOPT  = 3:  Prints  1 word  of  Hollerith  information. 

N should  equal  1 . 
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The  above  call  adds  the  N Hollerith  code  words  to  the  list  of  code 
words  to  be  output  on  the  next  line  of  print.  When  8 code  words  have  been 
accumulated,  the  line  is  printed;  any  excess  code  words  are  added  to 
the  list  for  the  next  line  of  print.  Codes  are  printed  with  "7X,A6,7(9X,A6)" 
format. 

(2)  The  printing  of  Values  of  Variables 

CALL  STOVAR  (N,  A1 , A?,  . . . , A8)  1 < N < 8 
where  N is  an  integer  identifying  the  number  of  arguments  following  it. 

If  N <_  8,  then  the  remaining  arguments  must  be  dummy  arguments.  For 
example,  CALL  STOVAR  (6,  A1 , A2,  A3,  A4,  A5,  A6,  DU,  DU). 

CALL  STFL  (JOPT,  N,  ARG1).  This  call  prints  the  value 
of  one  variable  when  JOPT  = 1,  N = 1,  and  ARG1  is  the  name  of  the 
variable. 

Lines  accounting  is  taken  care  of  within  this  routine. 

Values  are  printed  with  a "1PE15.7"  format. 

(3)  Forcing  final  print 

To  force  any  possible  remaining  print 
CALL  STFL  (0,  1,  DU) 
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FLOW  DIAGRAM  (DEF.) 


NPAGE  « NPAGE+I 
LONG  « t 
EJECT  PAGE 

write:  “six- degrees 

OF  FREEDOM  FLIGHT  PATH 
STUDY  GENERALIZED 
COMPUTER  PPOGRAM" 

WRITE  " INDSOF  Z" 
CASE  ( NCASE) ■ STAGS 
INSTAGE).  PAGE  (NPAGE) 
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' STFL  > 
, ENTRY  j 


J#PT 

I * O 


GETARG,  • ARG, 
WENT  * 5 
K • 2 

CALL  ARRAY <N,0) 


RETURN 


CLOW  DIAGRAM  (STFL) 


/STFLD\ 

[entry/ 


INTI ALIZES  BY 
DATA  STATEMENT: 
IZ  * O 

CLEAN  *.TRUE. 
INTEG  = FALSE. 


ARG2 

ARG  1 j 

...  N 

{ *91  > 

NENT  * I 


NENT  » 2 

/ N 

1 

' 

K * 1 

^return) 

CLEAN 

* TRUE.  1 

CALL  ARRAY (N.O) 

CALL  ARRAY  ( N,  1 > | 

(RETURN 


RETURN  Hs- 
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FLOW  DIAGRAM  (ST$VAR) 


/ST$VAR> 
I ENTRY  j 


I N | £ 2 
NO 

’ | 

I N | S 4 

I NO 


RETURN^ 
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10.  LINES  - LINES  ACCOUNTING  ROUTINE 

a.  Purpose  - To  keep  an  accounting  of  the  number  of  lines  printed 
per  page,  and  to  provide  for  page  control. 

b.  Method  - If  the  number  of  lines  to  be  printed  (LC0UNT)  is  such 
that  it  will  not  fit  on  the  current  pane,  the  page  is  ejected  (via  DEF) 
and  printing  will  begin  on  the  new  page.  Initially,  the  location  L0NG 
should  be  set  to  zero,  indicating  that  currently  no  lines  have  been 
printed  on  the  present  page. 

c.  Usage  - Entry  is  made  via  the  statement, 

CALL  LINES  (LC0UNT) 

where 

LC0UNT  = A fixed  point  variable  or  constant  indicating  the 
number  of  lines  to  be  printed. 

11.  ASIN  - ARC  SINE  FUNCTION 

a.  Purpose  - To  compute  the  arc  sine  of  a normalized  floating 
point  argument  (X). 

b.  Method  - , , 

Sin"1  (X)  = J - Cos-1  (X) 

c.  Usage  - The  Arc  Sine  is  computed  using  the  statement,  Y = ASIN  (X) 

where  jX|  < 1.,  and  Y = Sin"^  (X). 

12.  ACOS  - ARC  COSINE  FUNCTION 

a.  Purpose  - To  compute  the  arc  cosine  of  a normalized  floating  point 
argument  X. 

b.  Method  - For  |X|  < 7.4505806  x 10  the  arc  cosine  Is  set  equal 
to  tt/2.  For  X = 1.,  arc  cosine  is  set  equal  to  zero,  and  for  X = -1  arc 
cosine  is  set  equal  to  tt  . When  the  argument  X / + 1 . , the  routine  gives 
the  arc  cosine  in  radians  from  C to  tt. 
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FLOW  DIAGRAM  (ASIN.) 
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c.  Usage  - The  arc  cosine  is  computed  using  the  statement 
Y = AC0S  (X) 

where  |X|  £ 1.  and  Y ~ Cos"1  (X). 

13.  ATAN2  - ARC  TANGENT  ROUTINE 

a.  Purpose  - To  compute  the  arctangent  of  the  quotient  of  two 
normalized  floating  point  quantities  Y/X,  with  proper  quadrant  control. 

b.  Method  - The  routine  computes  the  quotient  Y/X.  The  arctangent 

is  computed  with  quadrant  according  to  the  sign  of  Y and  X.  If  X = 0 

-1  Y * 

and  Y f 0,  the  routine  computes  1 - Tan  (Y/X)  = j y j*  . 

If  X = 0 and  Y = C,  it  computes 
Z = Tan"1  (Y/X)  = 0. 

c.  Usage  - The  arctangent  of  Y/X  is  obtained  via 
Z=ATAN2(Y,X) 

14.  ERROR,  EXERR  - GENERAL  TABLE  ERROR  ROUTINE 

a.  Purpose  - To  provide  a method  of  indicating  the  table  which 
may  possibly  contain  an  error.  Also,  to  provide  the  Stop  Number  that 
identifies  an  error  condition  in  an  equation. 

b.  Method  - By  the  statement  CALL  ERROR  (LOCT)  in  which  LOCT  is 

the  subscript  of  the  curve  in  error,  the  routine  will  search  the  subscript 
table  and  find  the  corresponding  BCD  word.  This  word  will  then  be 
printed  as: 

"TABLE  ERROR  AAAAAA" 
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Where  AAAAAA  is  the  BCD  name  of  the  table.  If  the  name  cannot  be  found 
in  the  directory 

"TABLE  ERROR 

LOCATION  OF  TABLE  NOT  LISTED  IN  DIRECTORY"  is  printed  and  a 
return  to  the  calling  program  is  made.  In  either  case  INDSTE  is  set 
to  zero.  By  the  statement  CALL  EXERR  (NUM)  in  which  NUM  is  the  stop 
number,  the  routine  will  write  "STOP  NUMBER  III" 

where  III  = the  stop  number.  If  NUM  = 0,  an  exit  is  made  from  the  routine 
with  no  printing.  In  either  case  INDSTE  is  set  to  zero. 

c.  Usage  - Entries  are  made  to  the  routine  with  the  following 
statement: 

(1)  CALL  ERROR  (LOCT) 

where  LOCT  is  the  table  subscript. 

(2)  CALL  EXERR  (NUM) 

where  NUM  is  the  stop  number. 

15.  NDTLU  - N-DIMENSIONAL  TABLE  LOOK-UP  ROUTINE 

a.  Purpose  - To  provide  a method  of  linearly  interpolating  in  a 
table  of  n independent  variables. 

b.  Method  - Given  the  arguments  X(l),  X(2),  ....  X(N-l),  the 

routine  computes  Y = F(X(1),  X(2) X(N-l))  by  linear  interpolation 

from  a table. 


61 


AFFDL-TR-71-155 
PART  IV 


Pt-OW  DIAGRAM  (ERR^R,  EXERR) 
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c.  Usage  - Entry  is  made  via  the  statement: 

CALL  NDTLU  (ND,  NA,  X,  Z,  XA,  ZR,  IE) 

where 

ND  = Dimension  of  lock-up  (when  V = F(X),  ND  = 2) 

NA  - An  array  of  length  ND-1.  Numbers  of  values  of  each  table  of  X. 

The  tables  are  listed  by  size. 

X = Tables  of  each  X in  order. 


Z = Functi 

on  Vali 

ues. 

If 

array 

must  be  ir. 

the 

Assume  NX=4,  NY 

=3,  NZ; 

=2. 

W(l)  = 

F ( XI , 

Yl, 

Zl) 

W(2)  = 

F(X2, 

Yl, 

zi) 

W(3)  = 

F(X3, 

Yl, 

Zl) 

W(4)  = 

F(X4 , 

Yl, 

zi) 

W(5)  = 

F (XI , 

Y2, 

Zl) 

W(6)  = 

F(X2, 

Y2, 

zi) 

W(7)  = 

F(X3, 

Y2, 

Zl) 

W(8)  = 

F(X4 , 

Y2, 

Zl) 

W(9)  = 

F ( XI , 

Y3, 

Zl) 

W(10)= 

F(X2, 

Y3, 

Zl) 

W(ll)= 

F(X3, 

Y3, 

Zl) 

W(12)= 

F(X4, 

Y3, 

Zl) 

ZR  = Results 
IE  = Error  Code 
= 0 No  error 
-1  X array  toe  small 
1 X array  too  large 


F(X,  Y,  Z)  the  dependent  variable 
owing  order. 

W(13)=  F ( XI , Yl,  Z2) 

W(14)=  F(X2,  Yl,  12) 

W(15)=  F(X3,  Yl,  12) 

W ( 1 6 )=  F(X4,  Yl,  12) 

W(17)=  F(X1 , Y2 , 12) 

W(18)=  F(X2,  Y2,  Z2) 

W(19)=  F(X3,  Y2,  Z2) 

W ( 20 )=  F(X4,  Y2 , 12) 

W ( 21  )=  F ( XI , Y3,  Z2) 

W(22 )=  F(X2,  Y3,  12) 

W(23)=  F ( X 3 , Y3,  U) 

W(24)=  F(X4,  Y3,  7.2) 


2 array  not  ir.  ascending  order 
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Let  n be  number  of  indpendent  variables,  then  the  table  is  called  an  "(n+1) 
dimensional  table." 

Z = F(X1 , . . . , Xn) 

To  use  a table  of  dimension  1 3 ana  ± 5,  a call  to  HIH0  should  be 
made  with  the  list  of  arguments  in  the  calling  sequence  in  the  same  order 
as  the  independent  Variables  are  numbered. 

16.  ATMS  - ATMOSPHERE  CALCULATION  ROUTINE  (1969) 

a.  Purpose  - To  compute  the  atmosphere  characteristics:  Density, 
speed  of  sound,  pressure,  temperature,  and  kinematic  viscosity.  All 
are  a function  of  altitude. 

b.  Method  - All  atmosphere  characteristics  are  computed  using  the 
1969  ARDC  model  atmosphere.  Values  of  the  atmosphere  characteristics 
are  computed  for  positive  altitudes.  If  an  altitude  is  negative,  the 
sea  level  value  will  be  obtained. 

c.  Usage  - Linkage  is  affected  by 
CALL  ATMS  (HGC7F) 

where  HGC7F  = altitude  in  feet 

17.  INVR3  - INVERSE  OF  A NONSINGULAR  3X3  MATRIX 

a.  Purpose  - To  compute  the  inverse  of  a nonsingular  3x3  matrix. 

b.  Method  - Let  a = [a..]  i , j = 1,2,3 

• J 

then  A"'*  - (ii  [A..]  is  computed  where  A.,  is  the  cofactor  of  a... 
The  matrix  A must  be  stored  in  the  usual  FORTRAN  sense  ( i . e . , 
columnwise) . 
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FJLOWD I AG  RAM  (NDlLU) 
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FLOW  DIAGRAM  (NDTLU) 


RETUR 
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MB  ( 1 ) * 

0. 

MB ( 2 ) - 

11000. 

MB (3)  ■ 

25000 

HB(4)  = 

47000. 

MB  ( 5)  * 

53000. 

HB  ( 6)  = 

79000. 

HB  ( 7)  = 

90000. 

MB (8)  * 

105000. 

HB  ( 9 ) * 

160000. 

HB  (10)  = 

170000 

HB  (It)  ' 

200000. 

RH^B(l) 

= 23.7692E  - 4 

RH$B(2) 

= 7.0620  E -4 

RH0B ( 3 ) 

= 7.7650E  -5 

RHWB(4) 

= 2.8804E-6 

RH9B ( 5) 

= I3.9468E  - 7 

RH#B(6) 

= 4.II89E  -8 

RH$B(7) 

= 4.26I0E-9 

RH$B (8) 

= 2.2320E  -10 

RH^B (9 ) 

= 1 8450E-I2 

RH()B  (10) 

= 1 .3380E-I2 

RH(^B  (II) 

= 6.II30E-I3 

PB(1)  = 

2116.21695 

PB(2)  = 

472  .73 

PB ( 3)  = 

51.979 

PB(4)  = 

2.5155 

PB (5 ) * 

1.2181 

PB  (6)  = 

2.I08E  -2 

P8(7)  = 

2I.809E  -4 

PB ( 8)  * 

I5.562E  -5 

PB(9)  = 

75.578E-7 

PB(IO)  = 

58.954 E -7 

PB  (II)  - 

29.759E-7 

TB( 1 ) = 

518.688 

TB( 2)  = 

389.988 

TB(  3)  * 

389.988 

TB{  4)  = 

LOS. 788 

T8( 5)  = 

508.188 

TB(6)  = 

298.188 

TB(  7)  * 

298.188 

TB  ( 8)  * 

406,'I88 

TB(  9)  = 2386.188 

TB(IO)  =2566.188 

TB  (II)  =2836.188 

FLOW  DIAGRAM  (ATMS) 


AK1  (1)  = 

- 225569E-4 

AKi  (2)  = 

.0 

AKi  (3)  = 

.138466  E -4 

AKi  (4)  = 

.0 

AKI  (5)  = 

-.159202  E -4 

AKI  (6)  = 

.0 

AKi  (7)  = 

.24I458E  -4 

AKI  (6)  = 

. 886289 E -4 

AKI  ( 9)  = 

.T54341E -5 

AKI  (10)  = 

.350715  E -5 

AKI  (II)  = 

. 222I29E  -5 

AK  2 < 1 ) = 

- 5.25612 

AK  2 ( 2 ) = 

0. 

AK  2 ( 3 ) = 

11.3883 

AK  2 ( 4 ) = 

0. 

AK  2 ( 5 ) = 

- 7.59218 

AK  2(6 ) = 

0. 

AK  2 (7  ) = 

8.5412 

AK  2 ( 8)  = 

1.70824 

AK2(9) = 

3.41648 

AK  2 (10)  = 

6.83296 

AK  2 (1 1)  = 

9.76137 

AK  3 ( 1 ) = 

.0 

AK  3 ( 2)  = 

I57689E-3 

AK3  ( 3 ) = 

.0 

AK  3 (4 ) = 

.I20869E-  3 

AK3 (51 = 

.0 

AK  3(6)  = 

.206234E- 3 

AK3 (7) = 

.0 

AK3  (8)  = 

.0 

AK3  (9)  = 

.0 

AK3  (Id  = 

.0 

AK3  (II)  = 

.0 

All) 

.759511 

A (2 ) * 

.935787 

8(1)  = 

.174164 

B ( 2)  = 

.273966 

C(i)  = 

220. 

C (2 ) = 

180. 

D(i)  = 

25. 

0(2)  = 

140. 

r»v. 
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INDATM  * 0 


.MGC7F  > 2.5E6. 


FLOW  DIAGRAM  (ATMS) 


HGP  > 90000, 


VS77F  * 0. 
TA77R  = 0. 
PA77P  = 0. 
ANUA7F  = 0. 
RH$AS  = 0. 


RETURN 


HGT90  = . FALSE. 


HGP  < HB  (11) 


LAY  = 11 


HGC7F  < 0. 


VS77F  = 1116.43372 
TA77R  = TB (1) 

PA77P  = PB ( 1 ) 

ANUA7F  = I.3723288E-4 
RH^AS  = RH$B  ( 1 ) 


1=1  + 1 


,^i/H6P>HB  ( I ) 


LAY  = I - 1 


TMP  = .3048*  HGC7F 
HGP  = TMP/I.  + 
TMP/6356766.) 

HCT90  = TP'.'E 
m = a 


HGP  > 180000. 


i a m 


TMP  = HGP  - HB  (LAY) 
TMP2  = I.+  AK1  (LAY  * TMP 


AFFDL-TR-71 -155 
PART  IV 


FLOW  DIAGRAM  (ATMS) 


1 

TMP3  = EXP(  -AK3 (LAY) 
• TMP) 

PA77P  « PB(LAY)*  TMP3 
RM$AS  * RH<>B  (LAY)* 
TMP3 


1 


, 

F 

STA77R  = SQRT(TA77R) 
VS77F  = 49.020576 

* STA77R  # 

ANUA7F  = .02269886-6 
« (TA77R*  STA77R 
/((TA77R  + 198.72) 

• RH$AS)) 

VS77F  = 0. 

ANUA7F  = O. 

TA77R  - TA77R  * ( A ( M ) 

- B(M)  x ATAN  (KGP 

- C(M )/D(M ))) 


7 0 
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c.  Usage  - Linkage  is  obtained  via  the  statement: 

CALL  INVR3  (A,B,  INDER) 

where 

A = the  array  name  of  the  matrix  to  be  inverted. 

B = the  array  name  where  the  resulting  matrix  is  to  be  stored. 

INDER  = an  error  indicator  set  by  the  routine. 

(a)  INDER  = 1,  results  are  good 

(b)  INDER  = 2,  A was  0. 

18.  MULT31  - A MATRIX  MULTIPLICATION  ROUTINE 

a.  Purpose  - To  postmultiply  a 3 x 3 matrix  by  a 3 x 1 matrix. 

b.  Method  - The  result  of  [A]  [B]  = [C]  is  computed  using  single 
precision  floating  point  arithmetic.  All  elements  must  be  stored  in 
the  normal  FORTRAN  sense  (i.e.,  columnwise). 

c.  Usage  - The  matrix  multiplication  is  obtained  by  the  statement: 

CALL  MULT31 (A,B,C) 

where 

A = array  name  of  the  3x3  matrix  [A] 

B = array  name  of  the  3 x 1 matrix  [B] 

C = array  name  of  the  resulting  3 x l matrix  [C] 
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19.  TRNP0S  - A 3 x 3 MATRIX  TRANSPOSE  ROUTINE 

a.  Purpose  - To  transpose  a 3 x 3 matrix  [A]  to  obtain  he  3x3 
matrix  [A]' . 

b.  Method  - The  resulting  transposed  matrix  is  stored  in  a 
separate  array.  All  elements  of  [A]  must  be  stored  in  the  normal 
FORTRAN  sense  (i.e.,  columnwise). 

c.  Usage  - The  transpose  of  a 3 x 3 matrix  [A]  is  obtained  by: 

CALL  TRNP0S  (A,B) 

where 

A = The  array  name  of  the  3x3  matrix  [A]. 

B = The  array  name  of  the  3x3  matrix  [A]. 

20.  HIH0  - N-DIMENSIONAL  TABLE  CALL  ROUTINE 

a.  Purpose  - To  set  up  the  NA  array  and  Z location  of  tables  with 
dimension  from  3 to  5 as  required  by  the  calling  sequence  to  NDTLU, 
which  is 

CALL  NDTLU(ND,NA,X,A,XA,ZR,IE),  to  make  the  call  to  NDTLU,  and  return 
the  function  value  or  data  on  a table  read  error. 
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b.  Usaqe  - Linkage  to  the  subroutine  is  made  via  the  statement 
CALL  HIH0  (N,  L0CT,  NX1 , NX2,  MX 3,  NX4,  X1ARG,  X2ARG,  X3ARG,  X4ARG,  A) 
where 

N = dimension  of  table  look-up,  when  A = F(x),  N=2. 

L0CT  = location  of  the  first  value  in  the  table. 

NX1  to  NX4  - location  cf  number  of  points  in  the  XI  to  X4 
array  of  independent  variable  values 
X1ARG  to  X4ARG  - name  of  XI  to  X4  argument  or  a dummy  location  if  N < 5 
A = location  of  the  dependent  variable. 

21.  TLU  - TWO-DIMENSIONAL  TABLE  LOOK-UP  ROUTINE 

a.  Purpose  - Given  an  argument  X,  to  compute  Y = F(X)  from  a table 
of  X and  Y value:  <v  linear  interpolation. 

b.  Method  - The  table  of  X values  is  searched  until  for  some 
xi  < X < Xi+1 

Linear  interpolation  is  then  performed.  If,  for  some  i,  X=Xi  then  Y 
is  set  to  Y^ . 

c.  Usage  - Ent^y  is  made  via  the  statement, 

CALL  TLU  (X,  L0CT,  Y) 

where 

X = = variable  name  of  the  argument 

L0CT  = location  of  first  subscript  of  desired  table  of  data 
Y = variable  name  of  the  interpolated  value. 

All  tables  are  stored  in  a table  called  C.  Therefore,  i.he  first  value 
of  a particular  table  is  located  at  C(L0CT).  As  an  example,  assume 
foat  the  name  of  the  desired  table  is  FTAB01 . Then  the  curve  must 
be  stored  as  follows: 


C(L0CT)  = FTAB01  (1)  - N=iio.  of  points  In  curve  (integer) 
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( HIH(J  1 
iENTRY  , 


4 


XA  ( 1 ) 

= 

XI  AR6 

XA(  2) 

s 

X2  ARC 

XA(3) 

s 

X3  ARC 

XA<4) 

= 

X4  ARG 

NA  ( i ) 

NX  1 

NA(2) 

= 

NX  2 

NA'35 

= 

NX  3 

NA(4) 

= 

NX4 

NAT 

* 

NX  1 

IERR$R=  O 

Tno 


CALL^IL 
CALL  ERROR  (LfrCT) 
CALL  LINES  (2) 
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C(L0CT+1)  = FTAB01  (2)  = X, 
C(L0CT+2)  = FTABO'I  (3)  = V] 
C(L0CT+3)  = FTABOl  (4)  = X2 
C(L0CT+4)  = FTABOl  (5)  = Y? 


C(L0CT+2N)  = FTABOl  (2N)  = XN 
C(L0CT+2N+1 )=  FTABOl  (2N+1)  « Y^ 

22.  TFFS1  - ENGINE  THRUST  AND  THROTTLE  SETTING 

a.  Purpose  - To  provide  a method  of  introducing  the  engine 
thrust  characteristics  into  the  computation,  with  an  option  to  find  the 
throttle  setting  that  corresponds  with  a certain  thrust  and  Mach  number. 

b.  Usage  - Linkage  to  TFFS  is  accomplished  via  the  following 
statements: 

( 1 ) CALL  TFFS! 

Pre-data  initialization. 

(2)  CALL  TFFS3 

Thrust  computation  section. 

(3)  CALL  TFFS4 
Initial  print 

(4)  CALL  TFFS5 

Code  printing  to  identify  the  coming  time  history. 

(5)  CALL  TFFS6 

Time  History  Print 
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(6)  CALL  TFFS7 

Update  integration  (none  for  this  subprogram) 

23.  VPCS  - Vehicle  Physical  Characteristics 

a.  Purpose  - To  introduce  various  physical  characteristics  into 
the  general  solution  of  the  problem.  Here,  mass,  moments,  and  products 
of  inertia  of  the  vehicle  and  rotating  machinery,  reference  lengths  aid 
areas  for  aerodynamic  coefficients,  jet  damping  characteristics  lengtns, 
and  center  of  gravity  information  are  included. 

b.  Usage  - Linkage  to  VPCS  is  provided  ly  the  following  statements: 

(1)  CALL  VPCS1 

Pre-data  initialization.  Necessary  romina  values  are  set. 

(2)  CALL  VPCS2 

Post-data  initialization. 

(3)  CALL  VPCS3 

VPCS  computations  are  performed  if  INDVPC  is  non-zero. 
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FLOW  DIAGRAM  (TFFS1) 


SET  TO  zero:  T,  , 

Ty  , Tt  , Lf  , Mt  , Nt 
TXB7P  = 0. 
TYB7P  = 0. 
TZB7P  = 0. 
ALT77F  = 0. 
AMT77F=  0. 
ANT77F=  O. 


RETURN 


SET  TO  ZERO : T,  , 

Ty  | Tj  I Lj  , Mt,  N| 

TXB7P = 0. 
TYB7P = 0. 
TZB7P  * 0. 
ALT77F  * 0. 
AMT77F  = 0. 
ANT77F  * 0. 


^E2_(,nd7f^o") 

| N0 

I = 0 


1*1  + 1 

OBTAIN  T(I)-/(N(I),Mn): 
CALL  HIH<t> 

MT  (I)  ~ T(I)  » ZN  ( I ) 

NT  (I)  = -T(I)*  YN  (I) 
TXB7P  = TXB7P  + T (I) 
AMT77F  * AMT77F  + MT { I ) 
ANT77F  - ANT77F  +NT(  I ) 


j ' 

; j 


} 


i 

{ 

\ 
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(4)  CALL  VPCS4 
Initial  print 

(5)  CALL  VPCS5 

Code  printing  to  identify  the  coming  time  history  is  performed 
if  INDVPC  f 0. 

(6)  CALL  VFV6 

Time  history  print. 

24.  SACS  - AERODYNAMIC  FORCES  AND  MOMENTS 

a.  Purpose  - To  provide  a complete  accounting  of  the  various 
contributions  to  the  aerodynamic  forces  and  moments,  regardless  of  the 
flight  conditions  or  the  vehicle  considered.  In  SACS  coefficients 
are  computed  in  the  proper  coordinate  system  for  use  in  other  parts  of 
the  TOLA  program. 

b.  Usage  - Linkage  to  SACS  is  accomplished  via  the  following 
statements: 

(1)  CALL  SACS1 

Pre-data  initialization. 

(2)  CALL  SAC S3 
Aerodynamic  computation. 

(3)  CALL  SACS4 

Initial  print.  (None  for  this  subprogram). 

(4)  CALL  SACS5 

Code  print  to  identify  the  trajectory  point  for  SACS. 

(5)  CALL  SACS6 

Time  history  print  for  SACS. 

(6)  CALL  SACS 7 

Update  integration  (None  for  this  subprogram). 
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(7)  CALL  SACS8  (CLR,  a , CDR) 

Determine  the  a d that  corresponds  with  and  Cp. 

(8)  CALL  SACS 9 (o  , CL>  CD) 

Determine  the  C|_  and  CQ  that  corresponds  with  a . 

(9)  CALL  SACS10  (a  ) 

Determine  6 as  a function  of  a . and  all  current  variables, 
qri  d 

(10)  CALL  SACS11  (<$rN) 

Determine  6 „ as  a function  of  all  current  variables. 
rN 

25.  AER01  - AERODYNAMIC  DATA  LOOKUP  FUNCTION 

a.  Purpose  - To  look  up  aerodynamic  data  from  the  table  array  in 
C0MMON  TBDIR/C(300) . It  looks  up  the  first  two  values  of  a particular 
table  depending  on  the  argument  of  the  function. 

b.  Usage  - Linkage  to  AER01  is  accomplished  in  the  following  ways: 

Y = AER01  (L0CT,  AER02) 

where  L0CT  = subscript  for  the  table  C ( 300 ) 
which  will  be  the  1st  location  of  a particular  table  that  one  is  interested 
in. 

After  the  function  is  executed, 

AER01  = 1st  value  of  the  table  that  begins  at  C(L0CT) 

AER02  = 2nd  value  of  the  table  that  is  located  at  C(L0CT  + 1) 

26  . AQUAD  - A QUADRATIC  FUNCTION 

a.  Purpose  - To  solve  the  equation  a { x | x + b x + C = 0 

b.  Usage  - Linkage  to  AQUAD  is  accomplished  in  the  following  way: 

Y = AQUAD  (A,  B,  C,  XLLIM,  XULIM) 
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SACS! 

ENTRY 


FLOW  DIAGRAM  ( SACS  I ) 


WRITE  HEADING;" SAC r 


SET  7®  ZERO : a,nF , 
y ■ £ t w i n 

AA77P = o. 

AN77P  - 0. 

YA77P  = 0 
ALA77F  = 0. 
AMA77F  a O 
ANAZ7F  = 0. 


INDAa0=0  J 

NO 

WRITE  C®OE!  “CAVAH” 


WRITE  CDDES!  “CA 
CN  CY" 

WR.TE  COOES.'  "CL 
CM  CNN" 


~^RE1 

... 

YES 

1 NDAER=  0 

] 

ENTRY 

SACS4 


ENTRY 
SACS  7 


RETURN! 


® WRITE  VALUES  FOR!  CA, 
CN , CY 

• WRITE  VALUES  FDR;  CRM, 
CM, CYM 


WRITE  VALUE  FOR: CAVAH 


ENTRY 

SACSS 


| NO 

I NDAftO  = 0 


NO 

INDAER-0 
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YES 

CA  = CAMNU 
CN  = CNIMNU 
CY  s CYMNU 
CRM  =CLMNU 
CM  = CMMNU 
CYM  =CN2MNU 

0 

r~  — 

^ 1 

YA77P  : (EPSSn  CY 


+EPS6)  * DYNPP*  AREFF 
ALA77F  s (EPS7*CRM 
♦EPS )*  DYNPP 


* AREFF*DZRFF 
ANA?F  - ( EPSI  !*CYM 
+EPSI2)  * DYNPP 
•■AREFF  * D2RFF 
AA77P  = ( EPS3 • CA 
+EPS4) • DYNPP* AREFF 
ANA77P  * ( EPS  I * CN 
+EPS2)*  DYNPP*AREFF 
AMA77F r!EPS9»CM 
♦EPSlO)*  DYNPP 
• AREFF) • DIRFF 
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(SACS9) 

* 3 


0 


YES 


/ \ NO 

( I NDAO I = 0 )— *-  CA  = (AER®I(L©C(I),AER®2)-AER®2)»TMRI 
-| ' _ +AER®2 

c 


I NDA02 - 0 


YES 


C Aa*  = ( AERCI  ( LCC<2),  AERC2)  ~ AER®2)*Tw»P<l ) 

+AER®2 


INDAI 

YES 

l 

5=0  ^ 

1 

NO 

Cms  =(  AERCI  (L®C(I5),AER®2)-AERC3*TMP(I) 
9 +AER®2 

b 

1 

— 

[ N0  1 

1 1 

i 

YES 


/ * 

( I NDAI 6 = 0 V- •>  CNJs  2 = <AER®I(LCC(I6),  AER®2)-AER®2)*TMP\I) 

V S oq  +AERC2 


YES 


/ NO 

( I NDAIO-O  V- CN  = ' AERCI  (L©C(  10),  AERC2) -AER«2)*TMP(I) 
V ' +AERC2 


I NDAI  I =0  ^ 


NO 


YES 


c 


f 


CN  = (AERO I (L®C(  I I),  AER®2)  - AER®2)-*-TMP(l) 

+AERC2 


INDAi  2-0 


NO 


YES 


CN ~ .AERCI  (L®C(  12) , AER®2)  -AER®2  )*TMP(I) 
a +AER®2 


c 


NO 


INDA20  =0 


YES 


CN  * =(AER®i  (L®C(20),  AERC'2)  - AERC2)*  TMPO  j — 

+AERC2  1 


NO 

^ INDA22=0  J — ►- 


f" 

/ — 1 — \ r 

( INDA80  =0  ) — fc-J 


C.N  - ! AER®  I ( L®Cf  22) . AER®2) - AFR®2)»TMF<I ) 
**  +AERC2 


cAn  = (AERCI  ( L®C  (70) , AER®2)  -AERC2)»TMP(I) 

+AERC2 


0 


Y-^G)  (SA0S8I 
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INDA5I  - 0 
YEP  ~ 

INDA52  = 0 
YES  ~T~~  ~ 

_L LZ 

IN0A53=0 
YES  , 

INDA56=0 
YES  1 ^ ~ 


Cno  = (A£R®I(L®0(5I),AER®2/-AER®2)*TMP(I)+AER®2 
k = (AER© ! ( L«C  ( 52) , AER®2)  - AER®  2 ) • TMP(  I ) +AER®2 

C = ( AER®  I (UDC  (53) , AER®2)-AER{32  )*TMP(I)  +AER®2 

! — 

j C ,,gq  *(  AER©I(Lj®C(56)  , AER02)  -AER©2)#TMP(  I )+AER©2 


INDA57  = 0 )— Cmgq*=(AER®l(L®C(57),AER©2)-AERC2)*TMP<l)+AER®2 

I I 


INDA59  * 0 
~~YES  I 


I NDA58  = 0 J — »-j  Cma^-(AER®I(LCC(58),AER®2)- AER®2)»TMP(!)+AER®2 
YES  n * ’ 1 

\ NO 

( INDA59  * 0 V- *«-  Cmagg=(AER®l(ljDC(59.',  AE*®2)  - AER©2)ttTMP(l)+AER®2 

YES  | ~ *■  ~ — 1 

/ 1 NO  j 

( INDA60=0  Cm(ggq=(AER®Klj®C(60),AER®2)-AEW®2)*TMP(l)+AER®2 

YESC  1 ~ ~ 

^ * N NO 

f IMDA6I  = 0 >-*-  Cm(j  * ( AER®I  (LCC(6I),  AER3>2)-AER®2'*TMP(I)+AER®2 

YES  T ~ - 1 ~ — 

1 N0 

f INDA62  - 0 Cma  ( ” ( AER®  I ( L®C(62) , AER*>2)  - AER®2 ) *TMP(  I )+AER®2 


1 

f 

INDA62  - 0 

YES  1 

j 

INDA63- 0 ) — Cma  MAEROI  L©C(63),  AERfi?j- A£RC?!*rMP(l  ) + A£R©2 
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16  I (SACSIO) 


IN0A27  s C )“ H Cy£  - (A£R®!(L©C(27),AER®2)-AER«2)*TMP(I)+AER®2 


'Nf>A28=  0.  Cy^*  = (AER®l(L®C(28);AER(r2)-AER®2)»TMP(l)+AER®2 


INDA68-  0 Cnp  = {A£R©I(L©C(68),  AERC2)  - AER®2)*TMP;i)+AER®2 


INDA69-  0 V-^4  C,  0*  = (AER®I(L®C(69),AEP®2)-AER©2)*TMP(I)+AER®2 


INDA70-  0 V-*H  Cngr  - (A£R®I(L©C(70),  AER®2)-AER02;*TMP(I)+AER®2 


INDA7I  = 0 >-H  ^n8r2  s ( AER®  I ( L®C(  71 ) , AER®2 ) - AER©2)*TMP(  I)  +AER®2  \ 


17  J(SACSII) 


INDA03=  0 


CA£=(AER®I(L®C(3),  AER®2)-AER®2)*TMP(4)+AER©2 


IN0A04  =0  } — —4  CA^g*  = 'AER®I  (l®C(4), AER®£)- AER®2)*TMP(I)+AER®2 
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FLOW  DIAGRAM  (SACSl) 


, >- sNO 

7 INDA45  - 0 /-*- 


YES 

C INDA46 


0/ga  = ( AER® ! ( L®C  (45) , AER®2 ) - AER®2  )*TMP(I  )+A£R©2 


YES 


c 


N NO 

*°  f-*-  c/aSp" 


= ( AER® I ( L®C (46) , AER©2 ) - AER®2 )*TMP ( I ) +AER®2 


IN0A47  = 0 


NO 


YES 


C (L®C(4T)>AER®2J~AER®2)*TMP<I)+AER®2 


<*- 

INDA4 

1 N NO 

ie  = o *- 

C / = (AERO!  (L0>C(48),  AER<T2)-AER<D2)*TMP(I)+AER02 

— 

YES 

. NO  1 

i 

^ INDA49 - 0 y 


YES 


c 


C/f  = ( AEROI  (LfflC(49),  AER®2)  — AER®2  )-*TMP(l)  +AER02 


NO  f 


INDA50*  0 


YES 


K 


C / = ( AEROI  (L©C(SO),  AER®2)  -AER©2)*-TMP( I ) +AER®2 

r * 


h 


INDA54 


YES 


c 


N NO 

i °y~ 


Cm/g  = ( AF.R® I { L®C(54) , AER®2) -AER©2)*TMP( I ) +-AER®2 


NO 


INOA55  = 0 


YES 


C 


Cm(gZ  = (AF.R®l(L®C(55),  AER®2)-AER®2)*TMP(l)+AER®2j-  — . 


INDA64 


^rv: 


NO 


YES 


Cmqx  *<AER®Ul.®C(64),AER®2)-Afrt®2)*TMP(l>+AER©2 


^ INDA65=~Q  Y 

t ' 


Cno  = (AER©l(L®C(Si>),  ALR®2)-AER®2)«TMP(I)+AER®2 


^~INDA66^  0 ^ — »-|  Cnc 


: ( AER® I ( L©C l G6) , AER©2 ) - AER®2P* TMP(  I ) +AERf 
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Gy  = cyo  + cVa  lal  + cya2  a2  + Cy^/3  +cypz(3\(3\ 

♦Cygr  Sr  + Cygr2|Srl  Sr  + Cyagr|a|  Sr  + Ciafi\a\  (3 

+ cy£8r  1^1  Sr  + (cy/3  +cy/3xAxC.G.)/3d2/2VA 
+ (Cyr  + CyryAXQ  Q ) r d2  / 2 VA 

Crm  =G/0+  C*a  |ai  + C*a2  az  + Cip  |/3|  + C/gp  Sp 

+ci8p2  I^pI  Sp  + c*a8plal  Sp  + CW 3 la£| 

+ c/)SSpi/3l  SP  + C/pPd2/2vA  +(C/r+C/r)tAXc.G.)r02/2VA 

Cm  s Cm0  "*■  Cma  <J  + Gm<j2  a la!  + + @ 

+ CmSq  H CmSq2  Sp  I8pI  + Cm aSq  ial  8<* 

+Cm af3  a \/3\  + Cm££q  l#l  8p  + (Cmd  +Cma!£  AxC.G>  «<M2VA 

+(Cmq  + Cmq>i  AXCq  )qd,/2VA~CN  AXC  G /d. 


cym  = cr<o  < -r.a  lal  + cMa2  «a  + cn#£+<«vg£  + cn/98£l0l 
+cn8r  ^ + Cngr2!8r{  Sr  ^Cnagr!«ll  Sr  + Cna/Q|a|/3 

+Cn/38r  l#l  8r  ^ <cn/8 +co/8xAxC  G.)^d2/2VA 
+(Cnr  + Cnrx  AXC  G ) rd2  / 2VA— Cy  AXC  G./  d2 
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IAP  < 3 ] 

" ' 

NO 

|Hg-HR|<  10“*) 

' [no 


ENTRY 
SAC59  , 


IAP < 3 

| NO  ' 

|Hg-HR|<IO"# 

I NO 


1 G(D  1 =2 

— - 

— 1 

1 G<D  1 = 3 
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TMP2  » Cngr,  q*Sd2 
TMP3  • Cngr  q*Sd2 
TMP4  • [c„p-ry^  AXCG  /d2 

+ l/SdKCoyg*  -Cyp*  AXCG/d2l]  q*Sd2 


SOLVE  FOR  SrN  WHEKE  SrL<  SrN  S8ru: 
TMP2|8rN|8fN  + TMP38rN+TMP4*0 
(THIS  IS  DONE  BY  USING  THE  FUNCTION  AQUAD) 
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FLOW  DIAGRAM  ( AERO  I ) 


RETURN 
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where  A = the  coefficient  of  J X ) X 
B = the  coefficient  of  X 
C = the  constant  term 
XLLIM  = the  lower  limit  of  X 
XULIM  = the  upper  1 imit  of  X 

After  the  function  is  executived,  the  value  of  X will  be  stored  in  AQUAD 

27.  OPTI  - Six-Degree  -Of-Freedom  Trajectory  Program  Over  a Flat  Planet 
a.  Purpose  - This  program  permits  the  computation  of  the  six  rigid- 
body  degrees-of-freedom  of  a flight  vehicle  where  the  motion  is  assumed 
to  occur  over  a limited  portion  of  the  planet.  Under  this  assumption, 
the  motion  is  taken  to  occur  within  a rectangular  coordinate  system. 

This  program  of  computation  has  been  provided  to  aid  in  such  analyses  as: 

(1)  Boost-phase  dynamics. 

(2)  Pitch-roll -yaw  coupling  motion  investigations. 

(3)  Autopilot  response  to  parametric  disturbance. 

(4)  Landing,  approach,  and  flare-out  maneuvers,  etc. 

(5)  Landing  gear  dynamics. 

The  effects  upon  the  motion  of  the  body  of  the  following  items  may  also 
be  accounted  for: 

(1)  Three  components  of  wind  velocity  and  acceleration  measured 
in  a local-geocentric  Cartesian  coordinate  system. 

(2)  An  arbitrary  planetary  atmosphere. 

(3)  The  gyroscopic  torques  imposed  by  rotating  machinery  aboard 
the  vehicle. 
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(4)  The  action  of  an  arbitrary*  three-plane  autopilot  upon  the 
control  deflections  of  the  flight  vehicle  in  response  to  the  measured 
motion.  The  motion  may  be  obtained  from  a stabilized  platform.  In  this 
program  the  motion  is  more  conveniently  printed  out  in  the  Cartesian 
Coordinates, 
b.  Linkage 

(1)  CALL  0PT1 

Pre-data  initialization. 

(2)  CALL  0PT2 

Initialization  after  data  is  read  in  and  computation  for 
initial  time  point. 

(3)  CALL  0PT4 

Calculation  of  equation  of  motion,  and  all  other  variables 
that  are  desired. 

(4)  CALL  0PT6 

Time  history  printout. 

(5)  CALL  0PT7 

Update  variables  that  are  integrated  by  the  integration  routine. 

28.  LGEARi  - LANDING  GEAR  CALCULATIONS,  PART  I 

a.  Purpose  - The  LGEARI  subroutine  computes  the  effects  of  ground 
reaction  and  landing  gear  dynamics  on  the  motion  of  a landing  vehicle. 

A maximum  of  five  independent  landing  gears  may  be  used. 

b.  Usage  - Linkage  to  LGEARI  is  provided  by  the  following  statements: 
(1)  CALL  LGEARI 

Pre-data  initialization. 

* Currently,  only  a Pitch-Yaw-Roll  sequence  is  programmed,  due  to 
storage  limitations. 
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FLOW  DIAGRAM  ( $PTI  ) 


$FTi 

ENTRY 


INDSOF * 2 ] 

WRITE  COOES:  t,tf,  Xg, 
Yg,  h,u,y,  w,  p,  q,  r,  MN, 
VA,  q",  Xg,  Yg,  Zg 
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i 


k 


FLOW  DIAGRAM  (»TTI> 
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FLOW  DIAGRAM  (OPT!) 


ENTRY’ 
®PF4  . 


Zg  = -hg,  o"D  = °- 

ro-Y<  ai  s<r,y-,-r 


deg.t®  radians:  ,6,<p 


/,  = c«s  8 c®s  \i> 

4 - SIN  f 
t3  =-SIN  6 cos  f 

m,  * SIN0  SIN<£SIN^COS<£C®S0 
mz  s C®S^C®S<£ 

m,  * COS0  SIN<£+COS<jfrSlN0SIN<£ 

n,  * SIN0  Ct>S<p  + SIN  ^>SIN<^>COS0 
t.z  *-C®S<f/  SIN 

ns=  C®S0C®S<£-SIN<£SIN0SINt// 


Xg  = Vg  C®S  yxfiOS  Ofc 
Yg  = Vg  C®S  yQ  SIN  ofc 
Zg  =-Vg  SIN  /q 


X8S< 

N dt 

y«*J 

Pt  • 

Yg  dt 

0 

H 

' Zg  dt 

o /, 

J. 

*>|N 

v = m, 

m2 

m5  * 

Yg 

w n, 

nz 

ns 

Zg 

INDSKI 

P-  1 

1 



= ? S 9 REF 

* ms  gREF 

* n s 9 REF 

h = - 

Zg 

FORM  l ,m,n  MATRIX 
FROM  INTEGRATION 
TABLES 


*9  /.  U /s  u 

Yg  “ fH|  rrij  Wj  v 

2g  n,  nj,  ns  w 


Xg-j 

f’xg  dt 
0 

£ 

II 

pt  • 

Yq  dt 

j'zg  dt 

CALL  ATMS 
TC  COMPUTE 
P,  T.P.  Vbi,p 
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FLOW  DIAGRAM  (®PTI) 
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FLOW  DIA G RAM  ( 4> P T I ) 


ENTRY 

$PTC 


INOSKP ! 


PRINT  : t,  ttl  rig,  Yg,  h,  u,  v, 

P,  q.  r,  Mn,  Vq,  q*  Xg,  Yg,  Zq 


|-  --  {^hg>AMAXH  J 
,,  , YES 

I NOARC  * F INOARC -T 


I CALL  STGTST 


RETURN 


S>0>' 


INOAPC  * 0 

«■' 

YES  ^ 


PRINT:  a,  0 


INDGSR  * 0 

- ' 

YES  ^ 

INDFPA  = 0 
YES  ^ 

INOFPR  * 0^ 

YES  L 


PRINT  : Vn 


PR  I 


NT  . y,  <r 


PRINT  y,  cr 


r 

I 

INOSTE  * 0 

PRINT  : 0p,  ^p,<#.p  | 

T 

1 

3 

WRI 

'ST4>P- 


TE  : 
HGCTF" 


PRINT  : 
Aip,  Ayp,  Agp 


— — (lNOACM«Oy  - 2 PRINT;  ax,ay,ax 


:v 


RETURN 


i iNmt'CT  • nv_N® 

PRiPOT  ■ W*r 

PRINT  TIME  HISTORY 
FR0M  LGEAR.VPCS 
TFFS,  SACS 


INORMC  * 0^ 
I YES 


PRINT;  9, 


\ 
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FLOW  DIAGRAM  { <J) P T i > 
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(2)  CALL  LGEAR2 

Initialization  after  data  read  in. 

(3)  CALL  LGEAR3 

Landing  gear  dynamics  computation,  Part  I 

(4)  CALL  LGEAR4 

Initial  Print.  (None  in  this  subroutine). 

29.  LGEA3C  - LANDING  GEAR  CALCULATIONS,  PART  II 

a.  Purpose  - Continuation  of  the  calculations  of  the  effects  of 
ground  reaction  and  landing  gear  dynamics  on  the  motion  of  a landing 
vehicle. 

b.  Usage  - Linkage  to  LGEA3C  is  provided  by  the  following  statement: 
CALL  LGEA3C 


30-  SDFLGP  - PRINTING  OF  SDF  AND  LG  VARIABLES 


a.  Purpose  - The  SDFLGP  subroutine  saves  data  on  tape  for  plotting, 


prints  output,  and  updates  variables  for  the  integration  routine. 

b.  Usage  - Linkage  to  SDFLGP  is  provided  by  the  following  statements. 
(1)  CALL  SDFLGP 


Saves  data  on  tape  for  plotting. 


Input  cards  required: 


Coi umn: 


1 


o 

L 


IPLT  1 
I SDF  1 
ISTPL1  1 
ISTPL2  1 
IS1PL3  1 
ISTPL4  1 
ISTPL5  1 
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IPLT  = 1 denotes  that  data  will  be  saved  on  tape  unit  13  for  plotting. 
ISDF=!  denotes  that  rigid  body  data  will  be  saved  on  tape.ISTPL  denotes 
data  for  landing  gear  k will  be  saved  on  tape. 

(2)  CALL  LGEA6P 

Writes  code  on  tape  to  designate  last  record  for  a particular 
case  and  also  writes  EOF. 

(3)  CALL  LGEAR6 

Prints  output  of  the  Landing  Gear  calculations 

(4)  CALL  LGEAR7 

Updates  landing  gear  variables  that  are  integrated  by  the 
integration  routine. 

(5)  CALL  LGEAR5 

Prints  titles.  (Not  used,  only  a return  statement.) 
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FLOW  Qj  AG  RAM  (LGEARI ) 


HWMTM.r^T*^w"t"  w w? **?*’*?  rt»  *-*«"< 
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FLOW  DIAGRAM  (LGEARl) 


UO 
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FLOW  DIAGRAM  (LG 
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FLOW  DIAGRAM  (L  GEAR  I ) 


NSTRUT  NSTRUT  NSTRUT 

fTRX  1 Z FTRXj  , fTRY  ■ 1 ' FTRYj  • fTRZ  = Z FTR2j 

i = I i = I i = I 


ftra 

rL!I  rLI2  rL!3 

fTRX 

ftrb 

- 

rL2I  rL22  rL£3 

• 

ftry 

ftrc 

rL3I  rL32  rL33 

ftrz 

NSTRUT  NSTRUT  NSTRUT 

“TRX  * Z mTRM;  . “TRY  1 Z “TRY,  . “TR2  = Z “TRZ : 

i = I i = l i=l 


M.x 

rlii 

rLI2 

rLI3 

“TRX 

Mty 

rL2I 

rL22 

rL23 

■ 

“try 

Mth  ^ 

rL3I 

rL32 

rL33 

“trh 

L J 

NSTRUT 

fXM  = FTRA  + Z ">j  *i  o 1 31  . ^YM  s rTRB 
i - 1 


NSTRUT 

fZM  5 fTRC  + Z "i  »'■  °i  33 


i = l 
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© 
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FLOW  DIAGRAM  LGEA3C 
(i  = 1,2, NSTRUT  IN  Tills  SUBROUTINE  ) 


%i 

1 

°i  1 1 rLii  + °i'S  rL»  °ii  l 

RLl2  + °il3 

RLS2 

°ill  RLIS  + ail3  RLSS 

R>2,i 

Rl23i 

= 

RL  21 

RL22 

RL2S 

Rlsii 

Rl32j 

RlSSj 

°isi  RLll  + aiS3  RL»I  °isi 

RL,2  + °iS3 

RL32 

°ijl  RLI3  + ai33  RLS3 

RAX  j 

*1 

m, 

"i 

Rix 

+ Olj, 

( r Fj 

- Si  ) 

RAYj 

8 

/« 

n2 

Riy 

RAZj 

- 

4 

m. 

"S 

Ri* 

+ “'S3 

(rFj 

- Si  1 

^Oi 

= 

( Xg  - 

r8R  + 

> 

X 

RGS  + (Zg  + RAZj  ) 

RG,j 

Xr, 

- 

( Xg  - RgR  + 

RAXj)  RG| | + (Zg  + RAZj 

)RG 

13 

YRj 

r 

Yg  + 

RAYj 

FOR  i = 1,2, , NSTRUT 
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FLOW  DIAGRAM  LGEA3C 
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FLOW  DIAGRAM  LGEA3C 


Rt>*i  5 “ Si  a|33  "(rpi  “ q°l3|+  PR|y  “qRix 


RDXGj 

r • 1 

xg 

'l 

m,  n, 

RDX  j 

RDYGj 

s 

Yg 

+ 

f2 

mZ  n2 

S 

RD  Y j 

RDZG  j 

Zg 

'3 

m 3 n 3 

RDZj 

OH 

VTX  j = RG,|  RDXGj  + RG,3  RDZGj  -t-  WTj  RI22i  < r o j — S j ) 
VTYj  » ROYGj  - WT}  RI2,  j ( ro;  - 8; ) 

VTZj  * RG  3 j RDXGj  + RG33  RDZGj 
VGPTj  « ( V* TX j -*•  V* T Yj  ) i 
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FLOW  DIAGRAM  LGEA3C 


DXj  ' 

«G , , RAX; 

+ 

RGl3  RAZ, 

DY,  -- 

RAY; 

DZj  » 

F(G3I 

RAX 

i + 

Rt»33 

RAZj  + 

MTRX 

j * OY 

i FT 

RZ  i 

- OZj 

F T R ^ | 

MTRY 

i * DZ 

i FT 

R X j 

- OX; 

FT  RZ  : 
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FLOW  DIAGRAM.  LGEA3C 
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FLOW  DIAGRAM  (SDFLGP) 


SOFuGP ' 
ENTRY  , 


IPLT  = 0 
I NO 


(RETURN 


LASTPT  3 0 
IPPT  * IPPT  + I 
(SAVE  t IN 
AN  ARRAY) 


SAVE  THE  FOLLOWING 
VARIABLES  IN  ARRAYS: 

i Nm.  1,  ®p« 

'('p*  ^*P'  az> 

l>q,  Xq,  Yq,  a„,  Xq 


ISTPLI * 0 

~ ~™Tno 


SAVE  THE  FOLLOWING 
VARIABLES  IN  ARRAYS: 


PT I > SF|.  §|«  pl»  P2I' 
sl»  S|»  S|,  S2I>  S2I  » 

| s2|.  MAt.  *** T I . WTI 


ISTPL?.  * O 
YES 


SAVE  THE  FOLLOWING 
VARIABLES  IN  ARRAYS: 
pT2.SF2>  §2.  P2,  P22, 
S2>  Sg, Sg,  Sgg,  $22’ 
s22:  MA2 > **^T2 • w T2 


ISTPL3  = O 
I NO 


SAVE  THE  FOLLOWING 
VARIABLES  IN  ARRAYS: 
rT3.  Sp‘3. S3,p3.  p23, 
S3,  S3, S3,  S23.  S23, 
S23>  mA3>  “T3-  wT3 


ISTPL<  * O 
Tno 


SAVE  THE  FOLLOWING 
VARIABLES  IN  ARRAYS: 

PT4»  Sp4»  §4> p4»  p24> 
S4, S4, S4,  S24, S24, 
S24,MA4,u;T4,  wt4 


ISTPL5 * 0 
I NO 


SAVE  THE  FOLLOWING 
VARIABLES  ON  ARRAYS: 

PT5 > SF5-  ®5«  p5*  P,5, 
S5,  S5,S5,  S26,  Sg5« 

s25>  mA5»  ‘‘'TS-  wT5 
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FLOW  DIAGRAM  ISDFLfiPy 


33 


in  min 
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FLOW  DIAGRAM  (SOFLGP) 


WRITE  ARRAYS  OF  THE 
FOLLOWING  VARIABLES 
CN  TAPE  <IPPT  POINTS 
FOR  EACH  VARIABLE  I : 
fT5.  Sp5,  85, Pa.Pjj, 
S5, Sg, S5, S25,  ^25" 
S25,  MA5,  wr5,  wT5 


IPPT'O 


/TIME  LEFT 
ON  COMPUTER 
>20  SECONOS  ? 


YES 


NO 


LASTPT » I 

WRITE  LASTPT  ON  TAPE 
LTP  = 0 

WRITE  END-OF-PILE 
CN  TAPE  UNIT  13 


LTPT 1 0 
LASTPT  = I 

BACKSPACE  TAPE  UNIT 
13  ONE  RECORD 


WRITE  LASTPT  ON  TAPE 


rV/ — 

IRETURNH* LASTPT  * 


WRITE  END-OF-FILE 
ON  TAPE  UNIT  13 
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FLOW  DIAGRAM  ( SDR  LG  P ) 


ENTRY\ 

IlgearJ 


f ENTRY  \ 
Ilgear7; 


INOLG-0 


{RETURN) 


RETURN) 


INDLG'O 


WR!TE  :"LGEAR“ 


return) 


INOLG  *2 


INDLG  * 2 


WRITE  NAMES  AND 
VALUES  OF  THE 
FOLLOWING  VARIABLES: 
si.  Pi.  P2i.  PtI.Srj, 

[ sFi  i Oj.  Fczi.  /*i. 
VGPTi>  P'i'Rxi . FTRyi. 
[TRlj.  MAi,  M0j, 

88),  S, , Sj,  S,, 

^2i.  S2j,  S2j,  wTj, 
wTi  < • e I.  NSTRUT) 


write  names  and 

VALUES  OF  THE 
FOLLOWING  VARIABLES: 

FTi.  r*i.  FTR,i, 
FTRyi  - FTf}jr  j , 88  j , S j 
NSTRUT) 


UPDATE  THE  FOLLOWING 
VARIABLES  FOR  THE 
INTEGRATION  ROUTINE:  | 
$2i.  Szi 

• ' NSTRUT) 


UPDATE  THE  FOLLOWING 
VARIABLES  FOR  THE 
INTEGRATION  ROUTINE: 

^ i . S j , w-j-j 

■ I,  ....  NSTRUT) 


RETURN 


WRITE  NAMES  AND  i 

VALUES  OF  THE  ] 

FOLLOWING  VARIABLES: 
ftra.  fTrb.  fTrc.mt», 

MT,.*Titl  Fkw,  Fym, 

Fim * Lm> 


{RETURN 
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SECTION  V 
FORTRAN  EXTENDED 
OVERLAY  (1,0) 


1.  T0LAN1  - MAIN  PROGRAM  FOR  OVERLAY 

a.  Purpose  - A main  program  is  required  by  Fortran  Extended  for 

each  Overlay. 

b.  Usage  - Linkage  to  Overlay  (1,  0)  is  provided  by  the  following 
statement: 

CALL  OVERLAY  (T0LA1 , 1,0) 

c.  Flow  Chart 


call  diroda 
call  DiAiDA 
CALL  DIR2DA 


CALL 

READ 


CA 

TFf 

LL 

>S2 

1 1 
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2.  TFFS2  - SET  UP  TABLE  ROUTINE 

a.  Purpose  - To  call  TSRCH  to  set  up  table  for  various  subroutines. 

b.  Method  - When  this  routine  is  called  by  a certain  subroutine, 

it  will  in  turn  call  TSRCH  in  order  to  search  the  directory  for  the  exact 
BCD  name  required  by  the  particular  subroutine  that  called  this  routine. 

c.  Usage  - Entry  is  made  to  the  routine  with  the  following  statements: 

(1)  CALL  TFFS2 

Sets  up  tables  (or  provides  the  table  subscripts)  for  the 
subroutine  TFFS. 

(2)  CALL  SACS2 

Sets  up  tables  (or  provides  the  table  subscripts)  for  the 
subroutine  SACS. 

(3)  CALL  VPCTAB 

Sets  up  tables  (or  provides  the  table  subscripts)  for 
the  subroutine  VPCS2. 

(4)  CALL  OPTTAB 

Sets  up  tables  (or  provides  the  table  subscripts)  for 
the  subroutine  SDF2  (0PT2). 

(5)  CALL  LGTAB 

Sets  up  tables  (or  provides  the  table  subscripts)  for 
the  subroutine  LGEAR. 
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3.  READ  - INPUT  ROUTINE 

a.  Purpose  - To  provide  a general  method  of  reading  a variable 
field  data  card  and  assigning  variable  length  table.  Data  may  be  read 
into  symbolic  locations  in  memory. 

b.  Method  - Decimal,  Octal,  and  Integer  numbers  are  converted  to 
binary  integers.  BCD  information  is  stored  in  six-character  words. 

c.  Card  Formats  - The  variable  name  punched  in  Columns  1-6  is  the 
location  into  which  the  first  data  word  is  loaded.  The  variable  field 
information  is  for  relocation.  A fixed  point  integer  punched  anywhere 
in  the  field  (67-72)  will  be  treated  as  a subscript  to  the  variable  name 
punched  in  Column  1-5.  Negative  integers  punched  in  columns  67-72  will 

be  in  violation  of  the  subroutine.  The  first  blank  character  found  in  the 
card  to  the  right  of  Column  12  terminates  loading  from  che  card.  One 
exception  to  this  is  BCD  data.  Nine  6-character  words  may  be  loaded 
including  blanks. 

The  general  character  of  the  data  tc  be  loaded  is  determined  by  a 
three  letter  pseudo-operation  punched  in  columns  8-10.  The  pseudo- 
operations are:  DEC  or  blank,  0CT,  INT,  BCD,  and  TRA.  The  pseudo- 

operation TRA  is  a method  of  exit  from  the  subroutine. 

d.  Decimal  Data  - Decimal  data  beginning  in  column  12  and  ending  in 
column  66  is  converted  to  binary  and  loaded  into  the  symbolic  location 
punched  in  column  1-6  subscripted  by  the  integer  punched  in  column  67-72. 
Signs  are  indicated  by  + and  - preceding  the  number.  All  unsigned 
numbers  are  treated  as  positive.  If  either  the  character  E or  . or 
both  appear  in  the  decimal  data  word,  the  word  is  converted  to  a floating 
binary  number.  The  decimal  exponent  used  in  the  conversion  is  the 
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number  which  follows  immediately  after  the  character  E.  This  number  may 
have  a + or  - sign  preceding  it.  If  the  character  E does  not  appear 
the  exponent  is  assumed  to  be  2ero.  If  a decimal  point  does  not  appear 
it  is  assumed  to  be  at  the  right  of  the  numbers,  unless  it  is  the  only 
word  or  the  first  word  on  a card;  then  it  is  assumed  to  be  an  integer. 


All  the  examples  below  are  equivalent. 

(1)  12.345E03 

(2)  12.345E+03 

(3)  12.345E3 

(4)  1 2345E00 

(5)  12345. 

(6)  1.2345E4 

(7)  1 234500E-02 

(8)  +1 234500E-2 

Note  that  in  the  examples  above  all  decimal  words  have  decimal  points. 

If  the  frst  word  c.n  a card,  or  if  it  happens  to  be  the  only  word  on  a 
card,  does  not  contain  a decimal  point,  the  word  will  be  converted  to 
binary  integer. 

e.  Octal  Data-0ct  - The  Octal  data  is  loaded  the  same  as  decimal 

data  but  must  have  0CT  punched  in  columns  8,  9,  and  10.  All  data  is  converted 
to  binary  with  binary  point  assumed  at  the  right  end  of  each  word. 

f.  Hollerith  Data  - BCD  - Hollerith  information  is  loaded  from 
columns  3 through  66  and  assigned  consecutive  locations  for  every  6 
characters.  A maximum  of  nine  6-character  words  may  be  punched  on 
an.v  one  card  and  the  number  of  words  must  be  punched  in  column  12. 

A subscript  may  also  be  punched  in  columns  67-72. 
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g.  Transfer  - TRA--  The  purpose  of  the  TRA  card  is  to  transfer 
control  from  the  subroutine  back  to  the  main  program.  TRA  must  be  punched 
in  columns  8,  9,  and  10.  The  subscript  field  is  not  used.  A REWIND  may 
be  punched  beginning  in  column  12.  Only  the  R is  checked  and  the  only 
use  is  for  the  rewind  of  a data  tape. 

h.  Integer  - INT  - Integer  data  begins  in  column  12  and  ends  in 
column  66.  INT  is  punched  in  column  8,  9,  and  10.  It  may  be  relocated 
with  respect  to  the  BCD  name  by  punching  a subscript  integer  in 
columns  67  through  72.  If  only  one  data  word  is  punched  per  card, 
columns  8,  9,  and  10  may  be  left  blank. 

i.  Error  Messages  - A message  is  written  on  the  output  tape 
describing  the  type  of  error  encountered.  If  an  error  is  encountered, 
execution  of  the  case  is  deleted  and  the  subroutine  only  searches  for 
other  possible  errors  in  the  data.  The  following  error  messages  are 
possible. 

(1)  Symbol  not  in  directory. 

(2)  Column  12  is  blank. 

If  a bad  pseudo-operation  is  punched  in  columns  8,  9,  and  10  the 
subroutine  will  treat  it  as  decimal  data. 

All  checking  for  redundancies,  end  of  tape,  format  errors,  etc., 
ns  handled  by  FORTRAN  system  input/output  routines. 

j.  Usage  - Linkage  to  the  routine  is  made  by  the  following  statement: 
CALL  READ 
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k.  Data  Preparation  - The  first  card  expected  by  READ  is  a 
STCASE  TAB  with  the  S beginning  in  Column  1 and  TAB  punched  in  8,  9,  and  10. 
Following  this  card  is  a set  of  cards  which  define  the  table  sizes  necessary 
for  that  case. 

Example:  CTAB01  10 

CTAB02  20 

On  the  above  example,  CTAB01  is  punched  beginning  in  Column  1.  The 
numbers  10  and  20  are  punched  in  Column  12  and  indicate  the  number  of 
machine  cells  necessary  for  that  table.  Any  number  of  tables  may  be 
assigned  as  long  as  the  total  number  of  machine  cells  does  not  exceed  600. 
Follow  all  table  assignments  with  a TRA  punched  in  Columns  8,  9,  and  10. 

The  next  data  required  by  the  subroutine  is  a STCASE  and  followed  by 
any  combination  of  0CT,  BCD,  INT,  and  TRA  cards. 

4.  DSERCH  - DIRECTORY  SEARCH  ROUTINE  r0R  SUBSCRIPTS 

a.  Purpose  - To  provide  a method  of  searching  the  directory  to 
find  the  subscript  corresponding  to  a BCD  argument. 

b.  Method  - The  routine  searches  the  directory  for  the  exact  BCD 
name  required  by  the  argument.  When  an  equal  compare  has  been  found, 
the  corresponding  subscript  is  returned  as  a fixed  point  integer. 

c.  Usage  - Entry  is  made  to  the  routine  with  the  following 
statement: 

CALL  DSERCH  (SYM,  L0C , IER) 

where 

SYM  = BCD  name  being  search  for. 

L0C  = The  location  in  core 
IER  = Error  Code: 

IER  is  set  to  0 if  the  BCD  argument  does  compare  and 
is  not  a table  name. 
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IER  is  set.  to  -1  if  the  BCD  argument  does  compare  and 
is  a table  name. 

IER  is  set  to  +1  if  the  BCD  argument  does  not  compare. 

5.  TABRE  - TABLE  DIMENSION  SUBSCRIPT  ROUTINE 

a.  Purpose  - To  compute  subscripts  such  that  the  table  dimension 
requirements  may  be  variable. 

b.  Method  - Uses  input  data  prepared  by  the  user  to  compute 
subscripts  for  variable  table  assignments. 

c.  Usage  - Entry  is  made  via  the  statement 

CALL  TABRE 

When  the  READ  subroutine  processes  the  control  card:  STCASE  TAB 

then  the  subroutine  TABRE  is  called  by  the  READ  subroutine.  STCASE  begins 
in  Column  1 and  TAB  punched  in  Columns  8,  9,  and  10.  Following  this  card 
will  be  the  cards  requesting  table  sizes,  which  are  read  by  the  TABRE 
subroutine.  As  an  example  the  following  cards  may  be  read: 

TTAB10  10 

ATAB01  2 

TTAB10  and  ATAB01  punched  beginning  in  Column  1 and  the  required  machine 
cells  (10  and  2 in  this  case)  punched  beginning  in  Column  12.  Anything 
punched  past  Column  15  will  not  be  used.  After  all  table  assignments, 
a TRA  should  be  punched  in  Columns  8,  9,  and  10.  The  following  error 
messages  may  be  printed: 

1.  "Symbol  does  not  exist  in  table  list." 

2.  "Total  table  sizes  exceed  N,  change  maximum  total  table  sizes  to  NN 
in  subroutines  TLU,  HIH0,  TFFS,  AER0,  and  AUXR2,"  where  NN  is  the 
required  and  N is  the  maximum. 
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6.  TSRCH  - TABLE  SUBSCRIPT  SEARCH  ROUTINE 

a Purpose  - To  provide  a method  of  searching  the  directory  for 
table  subscripts. 

b.  Method  - The  routine  searches  the  directory  for  the  exact  BCD 

name  required  by  the  argument.  When  an  equal  has  been  found  the  corresponding 
subscript  is  returned  as  a fixed  point  integer. 

c.  Usage  - Entry  is  made  to  the  routine  with  the  following  statement: 
CALL  TSRCH  (SYM2,  L0C2,  N2,  IER) 

where 

SYM2  = BCD  name  of  argument 
L0C2  = Location  of  first  subscript 
N2  = NuPiber  of  sequential  subscripts  to  return  with 
IEK  = Error  Code: 

IER  is  set  to  +1  if  the  BCD  argument  does  not  compare 
IER  is  set  to  -1  if  the  BCD  argument  does  compare. 
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7.  ROUTINES  CALLED  BY  READ  SUBROUTINES 

The  following  are  subroutines  called  by  the  READ  subroutine  for  which 


flow  charts  are  not  provided. 
NAME  OF  SUBROUTINES 
DIPLAC 
READA 
STORE 
WRCARD 
PACKRR 
RITE 


LINKAGE 

CALL  DIPLAC  (RA1 , INC,  BLANK) 
CALL  PEADA  (IBCRW) 

CALL 'STORE  (N,  INC,  INX,  STAPE) 
CALL  WRCARD  (MSG) 

CALL  PACKRR  (II,  12,  NNN) 

CALL  RITE  ( IFI , FO,  JJ) 


8.  DlftOQA  - INPUT  DIRECTORY,  PART  I 

A routine  of  all  the  variables  in  data  statements  that  may  be  read 
by  the  input  routine  and  their  corresponding  subscript  lccation  in 
COMMON/ D I RCOM. 

9.  D I RIDA  - INPUT  DIRECTORY,  PART  II 
Continuation  of  DIRO. 

10.  DIR2DA  - INPUT  DIRECTORY,  PART  III 
Continuation  of  DIR1 . 
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SECTION  VI 
FORTRAN  EXTENDED 
OVERLAY  (2,  0) 

1.  T0LAN2  - MAIN  PROGRAM  FOR  OVERLAY 

a.  Purpose  - A main  program  ■'s  required  by  Fortran  Extended  for 
eacn  Overlay. 

b.  Usage  - Linkage  to  overlay  (2,  0)  is  provided  by  the  following 
statement: 

CALL  OVERLAY  (TOLA,  2,  0) 

c.  Flow  Chart 


147 


AFFDL-TR-71-155 
PART  IV 

2.  AUTS  - AUTOPILOT  AND  CONTROL  SYSTEM 

a.  Purpose  - The  AUTS  subprogram  calculates  the  control  responses 
that  are  needed  by  the  aerodynamic  forces  and  moments  subprogram. 

b.  Usage  - Linkage  to  AUTS  is  provided  by  the  following  statements: 
CALL  AUTS  - Compute  the  control  responses 

CALL  AUTSPR  - Compute  the  control  response  and  print  information 
as  indicated  by  input  parameters. 

3.  FLARE1  - AUTOPILOT  FLARE 

a.  Purpose  - The  FLARE  subroutine  calculates  the  actual  selected 
touchdown  conditions,  time  to  touchdown  for  X and  Y constraints,  required 
accelerations  to  meet  touchdown  conditions,  airspeed  flight  path  angle 
relative  tu  runway,  desired  thrust,  and  angle  of  attack. 

b.  Usage  - Linkage  to  FLARE  is  provided  by  the  following  statements: 
CALL  FLARE1  ( IPR)  Computes  all  that  is  mentioned  above. 

IPR  is  the  print  indicator.  If  IPR  = 0, 
the  print  routine  is  called, 

CALL  FLARE2  (IPR)  Only  computes  airspeed  flight  oath  angle 

relative  to  runway,  desired  Euler  roll  angle 
and  desired  angle  of  attack. 
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FLOW  DIAGRAM  (AUTS.i 
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4.  AUTPR1-  AUTOPILOT  PRINT  ROUTINE 


a.  Purpose  - The  AUTPR  subroutine  prints  output  for  the  autopilot. 

b.  Usage  - Linkage  to  AUTPR  is  provided  by  the  following  statements: 

(1)  CALL  AUTPR1  Write  !AUTS". 


(2)  CALL  AUTPR2 

(3)  CALL  AUTPR3 

(4)  CALL  PR4 

(5)  CALL  AUTPR5 

(6)  CALL  AUTPR6 

(7)  CALL  AUTPR7 


If  AUXICP  f 0,  WRITE 

"AUXILIARY  COMPUTATIONS" 

and  the  following  variables:  XR,  YR,  ZR, 

XR,  Yr,  Zr,  <|y 

If  MANL0G  f 0,  write  "GLIDE  SLOPE"  and  the 

following  variables:  Ve,  Vad,  yj,  QR,  CL  , 

R 

ad’  LR’  DR’  TTd’  HGS*  RR*  hea’  he*  heT5  hPa’ 
hPT*  ♦d 

If  MANL0G  f 0,  write  "FLARE",  and  the  following 
variables : t^,  Xjq,  ^TD* 

a • 

XRF.  Xrf,  hRF,  hRF,  Axr,  AhR,  ad,  TJd , <J>fl 
If  MANL0G  t 0,  write  “LANDING  ROLL",  and  the 
following  variables:  ISS,  ILR,  ICS,  BS,  Tj,  tr 

if  PITCHP  f 0,  write  "PITCH  AUTO  PILOT,"  and  the 
following  variables:  6 c ,,  a,  ad, 

ot  t » <5 

eT’  qd 

If  YAWAUP  f 0,  write  "YAW  AUT0PIL0T" , and  the 
following  variables;  SrN>  6,  6,  BeJ,  6rd>  U>et 

If  R0LLAP  f 0,  write  "ROLL  AUT0PIL0T" 
and  the  following  variables: 

^e*  ^eT’ 
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(8)  CALL  AUTPR8 


(9)  CALL  AUTPR9 


(10)  CALL  AUTP10 


(11)  CALL  AUTP11 


If  THR0AP  f 0,  write  "THR0TTLE  AUT0PIL0T" , and 

the  following  variables  N.. , I.  , i = 1,  2, 

a ai 

....  IN 

If  BRAKAP  f 0,  write  "BRAKE  AUT0PIL0T"  and  the 
following  variables:  , i = 1,  2, 

NSTRUT.  Also  if  BRAKAP  f 0 and  IBS  = 1 , 
write  the  following  brake  autopilot  variables: 

utri»  ^Ri*  wi£i » 1 = 1 » 2»  NSTRUT 

If  C0NTRP  t 0,  write  "C0NTR0L  RESP0NSE" , 

and  the  following  variables:  6„,  5 , 6 , 

q r p 

6q»  6r,  6p,  Ni . Ki , i = 1 , 2 IN. 

Also,  if  I ND I CP  f 0,  write  "INDICAT0RS" 
and  the  following  variables:  IAP,  ICi , 
i = 1,  2,  3,  4;  IBi#  i = k,  2,  3,  4,  5. 

If  ‘1ANL0G  f 0,  write  "FLARE" , and  the 
following  variables:  ,,d>  0rf,  TJd 


5.  THAUTS  - AUTOPILOT  THROTTLE 

a.  Purpose  - To  compute  the  desired  setting  and  desired  thrust  for 
each  engine. 

b.  Usage  - Linkage  to  the  throttle  autopilot  is  made  by  the  following 
statement. 

CALL  THAUTS  (I PR) 

If  IPR  = 0,  the  print  routine  is  called  to  print  the  throttle 
autopilot  output  variables. 
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NO 
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*R’  °R*  CLR*  ad' 

Lr,  Dr,  TYd,  h6s, 

**R»  h*a>  h*T* 
hPo • hPT-  4> d 


^MANlIg  * <TV^>-/rETURn) 

— n^r 


WRITE  : FLARE 

WRITE  FLARE  I 
VARIABLES  : »*,  1h, 

xTO’ hTD>  vxTD- vhTD> 
*RF>  *RF • hRF-  **RF * 
axR.  AhR.  <*d.  TTd. 

4>d 
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FLOW  DIAGRAM  (AUTPRI) 
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6.  ENGREV  - AUTOPILOT  ENGINE  REVERSE  LOGIC 

a.  Purpose  - ENGREV  is  a function  that  determines  if  an  engine 
is  to  be  reversed. 

b.  Usage  - Linkage  to  ENGREV  is  provided  by  the  following  statement 
Y = ENGREV  ( IR,  N,  Ng.  N,_R) 

where 

IR  = engine  reverse  capability 
N = engine  throttle  setting 
N„  = engine  throttle  constraint 

D 

^LR  = reverse  throttle  setting  for  reverse  signal  on  landing 

7.  ENGFL  - AUTOPILOT  ENGINE  FAILURE  LOGIC 

a.  Purpose  - The  ENGFL  subroutine  contains  the  logic  to  determine 
if  an  engine  has  failed  or  not. 

b.  Usage  - Linkage  to  ENGFL  is  provided  by  the  following  statement: 
CALL  ENGFL 

8.  CTENGL  - AUTOPILOT  COMMON  TWO-ENGINE  LOGIC 

a.  Purpose  - The  CTENGL  routine  calculates  the  distribution  of  a 
totei  dosireci  thrust  load,  T , between  two  engines,  A and  B,  as  T.. 
and  Tdg  respectively.  It  al.o  determines  the  corresponding  desired 
throttle  settings,  NdA  and  Ndg,  by  calling  CENGL.  twice. 

b.  Usage  - Linkage  to  CTENGL  is  provided  by  the  following  statement: 
CALL  CTENGL 
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9.  CENGL  - AUTOPILOT  COMMON  ENGINE  LOGIC 

a.  Purpose  - The  CENGL  routine  calculates  the  variable  N^  which 
is  required  by  throttle  autopilot  (THAUTS) 

b.  Usuge  - Linkage  to  CENGL  is  provided  by  the  following  statement: 
CALL  CENGL 

10.  TFFS8  - THROTTLE  SETTING  SEARCH  ROUTINE 

a.  Purpose  - Search  through  ranges  of  throttle  settings  to  find  the 
throttle  setting  TN  that  corresponds  with  current  thrust  TC  and  Mach 
number. 

b.  Usage  - CALL  TFFS8  (TC,  TN) 

11.  TFFS9  - COMPUTE  THRUST  as  F (MN,  TN) 

a.  Purpose  - Compute  thrust  as  a function  of  current  Mach  number 
and  designated  throttle  setting. 

b.  Usage  - CALL  TFFS9  (THRSET,  THRUST) 
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SECTION  VII 

PLOT  TAPE  GENERATING  PROGRAM  (PLTSDF) 

a.  Purpose  - To  generate  a plot  tape  from  a data  tape  that  was 
generated  by  the  Takeoff  and  Landing  Analysis  Computer  Program.  The  plot 
tape  is  plotted  by  the  CALCOMP  plotter. 

b.  Usage  - After  a plot  tape  has  been  generated  by  the  TOLA  program 
(see  subroutine  SDFLGP  for  data  setup),  this  program  is  turned  in  as  a 

s pa rate  job. 

c.  Input:  (1)  TAPE  - Data  generated  by  TOLA  program. 

(2)  CARDS  - The  following  may  be  read  in  as  input  using 
the  NAMELIST  feature  of  F0RTRAN  EXT.  If  the  value  of  any  variable  is  the 
same  as  its  nominal  value,  it  is  not  necessary  to  read  it. 


Name  Definition  Nominal  Value 

NCASES  = Number  of  sets  of  data  or  cases  to  be  plotted  1 

ISDFR  = 1 Rigid  body  data  is  stored  on  tape  1 

= 0 Rigid  body  data  is  not  stored  on  tape 
ISDF  = 0 Do  not  plot  rigid  body  data  0 

= 1 Plot  rigid  body  data 

ISTPR(l)  = 1 Landing  gear  data  for  gear  1 stored  on  tape  1 


= 0 Landing  gear  data  for  gear  1 is  not  stored  on  tape 
ISTPL(I)  = 0 Do  not  plot  landing  gear  data  for  landing  gear  1 0 

= 1 Do  not  plot  landing  gear  data  for  landing  gear  1 
IL  = 1 Do  not  plot  lower  chamber  pressure  and  2nd  Piston  0 
plots  (2nd  piston  was  not  uspd) 

= 0 Do  plot  lower  chamber  pressure  and  2nd  piston  plots 
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Name 

TFIRST  = 
TLAST  = 


PLTINT  = K 
FCTR  = 


XL 

YL 


Definition  Nominal  Value 

Trajectory  time  to  begin  plotting  0. 

Trajectory  time  to  stop  plotting  (Note=  0. 

if  both  TFIRST  = TLAST  = 0.,  the  entire 
time  history  on  tape  will  be  plotted) 

Plot  every  point  1 

The  current  factor  all  coordinates  are 
multiplied  by.  (That  is,  the  plot  is  made 
larger  or  smaller  if  FCTR  is  greater  than  1. 
or  less  than  1.  For  example,  if  wish  plot 
to  be  25%  of  the  original  size,  let  FCTR=.25) 

Length  of  X - Axis  in  inches  7.2 

Length  of  Y - Axis  in  inches  5.0 


Examples  of  input  data  cards 

Example  No.  1:  Plot  rigid  body  variables  and  landing  gear  variables 

for  gears  1,  3,  5.  Plot  every  point  and  plot  entire  time  history.  Assume 
rigid  body  and  landing  gear  variables  for  gears  1,  3,  5 are  stored  on  tape. 
EOR  (7/8/9) 


SINPUT  ISDF=1 , ISTPR=1 ,0,1 ,0,1 ,ISTPL=1 ,0,1 ,0,1$ 

E0J  ( 6/ 7/8/ 9) 

Example  No.  2:  Plot  landing  gear  variables  for  gear  No.  3.  Plot 
every  other  point  from  time  = 4.  to  time  * 10.  seconds.  Assume  rigid 
body  and  landing  gear  variables  for  gears  1,  2,  3,  4>  5 are  stored  on  tape 
for  time  = 0.  to  20.  seconds. 

EOR  (7/8/9) 

$ INPUT , ISTPL(T)=1 ,TFIRST=4. ,TLAST=10. ,PLTINT=2$ 

E0J  (6/ 7/8/ 9) 
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Example  No.  3:  Plot  rigid  body  variables  and  landing  gear 
variables  for  gear  No.  5.  Plot  every  point  and  plot  entire  time  history. 
Assume  rigid  body  and  landing  gear  variables  for  gears  1,  3,  5 are  stored 
on  tape,  with  the  size  of  graphs  to  be  50%  of  the  original  size  where  the 
original  size  of  the  X-Axis  is  8 inches  and  the  Y-Axis  is  6 inches. 

EOR  (7/8/9) 

$ INPUT , ISDF=1 , ISTPR-1 ,0,1 ,0,1 ,ISTPL(5)=1 ,XL=8.,YL=6. ,FCTR=.50$ 

EOJ  (6/7/S/9) 


d.  Output 

TAPE.  A plot  tape  is  to  be  plotted  by  the  CALC0MP  plotter.  The 
plot  tape  may  be  generated  to  give  the  following  plots: 

(a)  Rigid  body  plots 

Time  (sec)  Vs  Ground  Reaction  roll  moment  (lb>^0 
Vs  Ground  reaction  Ditch  moment  (Ib-ft) 

Vs  " " yaw  moment  (Ib-ft) 

Vs  Pitch  rate  (rad/sec) 

Vs  Yaw  angle  (Deg) 

" Vs  Roll  angle  (Deg) 

Vs  Z-Axis  ground  reaction  (lb) 

Vs  Z-Axis  acceleration  (ft/sec^) 

Vs  Mass  Center  sink  rate  (ft/sec) 

Vs  Mass  center  altitude  (ft) 

Vs  Cross  range  (ft) 

Vs  Mass  Center  altitude  (ft) 

Vs  Cross  range  (ft) 

Vs  X-Axis  acceleration  (ft/sec^) 

Vs  X-Axis  Velocity  (ft/sec) 


Down range  (ft) 

ii 

Time  (sec) 
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(b)  Landing  gear  plots  for  each  gear  (1) 

Time  (sec)  Vs  ground  reaction  along  strut  for  gear  No,  1 (lb) 
Vs  strut  resistance  for  gear  No.  1 (lb) 

" Vs  tire  deflection  for  gear  No.  1 (ft) 

2 

" Vs  upper  chamber  pressure  for  gear  No.  1 (lb/ft.  ) 

" Vs  lower  chamber  pressure  for  gear  No.  1 (lb/ft.^) 

" Vs  strut  acceleration  for  gear  No.  1 (ft/sec  ) 

" Vs  strut  velocity  for  gear  No.  1 (ft/sec) 

" Vs  strut  displacement  for  gear  No.  1 (ft) 

" Vs  2nd  piston  acceleration  for  gear  No.  1 (ft/sec  ) 

" Vs  2nd  piston  velocity  for  gear  No.  1 (ft/sec) 

" Vs  2nd  piston  displacement  for  gear  No.  1 (ft) 

Vs  wheel  axle  moment  for  gear  No.  1 (lb/ft) 

2 

Vs  tire  angular  acceleration  for  gear  1 (rad/s r „ ) 

" Vs  tire  angular  rate  for  gear  No.  1 (rad/spc) 

(c)  The  following  routines  are  calle*.  by  SDFPLi: 


PL0T 

SCALE 

PL0TE 

PL0TS 

LINE 

READ 

FACTOR 

SYMB0L 

BKSFIL 

AXIS 

(d)  The  appropriate  control  cards  are  required  after  the  job  card 

to  mount  tape  xxxx  from  TOLA  on  tape  unit  13  and  to  mount  tape  yyyy  for 

FLTSDF  output  (556  BPI)  on  tape  unit  7. 

(e)  If  TOLA  and  PLTSDF  are  run  on  the  same  job,  it  ’ c not  necessary 

to  use  tape  unit  13  for  storing  data.  In  this  case  disk  storage  may  be 

used. 
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